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#an se propose dans cet article de décrire quelques propriétés élémentaires 
(®, des équations aux dérivées partielles (e.d.p.) linéaires du second ordre à 
coefficients constants, autrement dit, dans le cas de deux variables, des équa- 
tions de la forme : 


d2u d2u d2u au au 
da Pay tua donhaptie EE à IE 


où a, b, c, a, B, y désignent six nombres réels donnés (a, b, c étant non tous nuls), 
F une fonction continue de deux variables réelles définie sur un ouvert U du plan 
et u une fonction inconnue, supposée de classe C 


On distingue a priori deux types de problèmes : 
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— ceux dans lesquels n'intervient pas la variable temps t, et qui ne dépendent 
donc que des variables spatiales x, y, z; ils sont appelés problèmes 
stationnaires ; | 


— ceux dans lesquels intervient, en plus des variables spatiales x, y, z, la varia- 
ble temps t ; ils sont appelés problèmes d'évolution. 


Æ On recherche le plus souvent des solutions vérifiant des conditions aux limi- 
tes, signifiant que la solution considérée u, a priori définie sur l’ouvert U du plan, 
satisfait certaines conditions sur la frontière de U. On distingue à ce sujet deux 
types de conditions, celles de Dirichlet et de Neumann. 


Les conditions de Dirichlet imposent à la solution u d’être continue sur l’adhé- 
rence de U, c’est-à-dire sur U et sa frontière, et d'être alors égale à une fonction 
donnée sur la frontière de U. 


Les conditions de Neumann imposent à la solution u d'être continue sur 
l'adhérence de U, c'est-à-dire sur U et sa frontière, et d'admettre en tout point de 
la frontière de U une dérivée du/aoN suivant le vecteur normal N orienté vers 
l'extérieur de la frontière de U (supposée suffisamment régulière) égale à une 
fonction donnée. : 


Dans un problème d'évolution, on recherche de plus des solutions vérifiant 
certaines conditions initiales (ou conditions de Cauchy), signifiant que, à l’ins- 
tant t = 0, la solution u(x, y, z, t ) de l'équation vérifie 

u (x, y, z, 0) = (x, y, 2) 


où f est une fonction donnée, et parfois 


9 
5 (X y, Z, 0) = g(x, y, x) 


où g est une fonction donnée. 


& Les problèmes que l’on peut alors étudier sont les suivants. 


e Un problème stationnaire donné avec des conditions aux limites ou un pro- 
blème d'évolution donné avec des conditions aux limites et des conditions initia- 
les, admettent-ils une solution et une seule ? 


e Dans l'affirmative, la solution obtenue dépend-elle continüment des don- 
nées (autrement dit, une « petite » erreur commise sur les conditions aux limites 
ou sur les conditions initiales conduit-elle à une «petite » erreur sur la 
solution) ? 


Notons dès maintenant que la linéarité de l'équation E implique que : 


— les solutions de l'équation homogène (équation obtenue lorsque le second 
membre F est nul) forment un espace vectoriel ; 


— la solution générale de l'équation complète s'obtient comme la somme 
d'une solution particulière de l'équation avec second membre et de la solution 
générale de l'équation homogène. 


Cet article est introductif et ne fait appel qu'à des méthodes élémentaires. En 
particulier, il ne fait jamais référence à la théorie des distributions, cependant 
centrale dans toutes ces questions. De même, les méthodes numériques de 
résolution (par différences finies ou par éléments finis) dont l'importance est 
essentielle puisque la résolution analytique n'est pas praticable en général ne 
sont pas abordées ici. 


On se reportera en bibliographie aux références [1] à [5]. 
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1. Classification des e.dl.p. 
linéaires du second ordre 


Ce paragraphe est destiné à distinguer trois types d'équations, 
qui se révèlent différentes tant du point de vue mathématique 
(propriétés des solutions, méthodes de démonstration) que 
physique. 

Étudions tout d’abord le cas des e.d.p. dépendant de deux varia- 
bles réelles. 


Définition 1. 
L'équation aux dérivées partielles (E) donnée dans l'intro- 
duction : 
CET dét a QU OU, gOU Ne Et y (1) 
dx2 oxdy ay? ox dy 
est dite de type : 
— hyperbolique lorsque A= bo 4ac>0; 
— parabolique lorsque A= b° —-4ac=0; 
— elliptique lorsque A=b?-4ac<0. 
Dans la suite, on dira que À = b? - 4ac est le discriminant de 
l'équation (1). : 


Cette définition est intéressante car invariante par changement de 
bases dans le plan, comme on le vérifie en effectuant le changement 
de bases, défini par : 


x'= Ax+ By ; 
y'=Cx+Dy ; 
avec AD-BCz0, 
et en posant alors: 
u'(x',y')=u"(Ax+By , Cx+ Dy) = u (x, y). 


En calculant les dérivées partielles de u en fonction de celles de 
u’, on voit alors que u est solution de (1) si, et seulement si, u” est 
solution de : 


; qu’, ; Du 2 du, à au”, g' 2u 
ax’? ax'à y’ dy? ox’ y” 


a 


7 (x, Y'X2) 


où l'on a notamment : 
a' = a? + bAB+ cB?; 
b' = 2(aAC + cBD) + b(AD + BC) ; 
c' = aC? + bCD + cD?. 


On vérifie par le calcul que le discriminant de (2) est égal à 
(AD - BC }? (b? - 4ac). Il est donc du signe de be - 4ac et les équa- 
tions (1) et (2) sont bien de même type. 

Proposition 1. 


A l’aide de changements de bases convenables et quitte à 
changer les notations, la recherche des solutions de l'équation aux 
dérivées partielles (1) est équivalente à la recherche : 


1. Des solutions de l'équation suivante dans le cas hyperbolique : 


o2u o2u 
de de TE Lequ= FN) 
ou 
o2u 
x 0y" * Bague FO 
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2. Des solutions de l‘équation suivante dans le cas parabolique : 


o2u u Lou v 
El: M LE Us F(X y) ; 
(On exclura le cas B = 0, qui ramène à une simple équation diffé- 
rentielle). 
3. Des solutions de l’équation suivante dans le cas elliptique : 
d2u o2u 
Er Er EL Er + yu = F (X y): 


Démonstration. © Établissons, par exemple, ce résultat dans le 
cas d’une équation elliptique. 

D'après les calculs de changements de base effectués précédem- 
ment, il suffit de trouver 4, B, C, D'tels que : 


a? + bAB + cB? = aC? + bCD + cD? ; 
2(aAC + cBD) + b(AD + BC)=0 
La première de ces relations est vérifiée en choisissant : 
°B=D=1; 
b 


° À, Cde la forme ph et —-h. 
2a 2a 


La seconde équation est alors vérifiée pour h? = (4ac - b?}/4&, 
expression positive puisque À < O ici. 

Comme AD-BC=2hZ0, ce changement de variables 
ramène clairement une équation elliptique à la forme indiquée 
dans la proposition 1. (0) 

Proposition 2. 

A l’aide de changements de fonctions inconnues convenables, la 
recherche des solutions de l'équation aux dérivées partielles (1) est 
équivalente à la recherche: 


— des solutions de l'équation suivante dans le cas hyperbolique : 


d2u d2v 


È- op + kv= G(X, y) 
ou ed: ÿ) ; 
OX dy 


— des solutions de l'équation suivante dans le cas parabolique : 


La k 2V 


ep G(x, y) ; 


— des solutions de l'équation suivante dans le cas elliptique : 


dv, d2v 


= Mer à +kv = G(X,y) 


Démonstration. © Reprenons encore le cas d'une équation 
elliptique qui, d'après la proposition 1, peut se ramener à la forme : 


22 d2u 
— ST 48 S J'u= FN. 
Posons alors 


U(X, y) = V(X y) exp(- BC EY), 


Un simple calcul de dérivées montre que l'équation précédente 
équivaut à la suivante :: 


dv, dv : ax + By 
Sat 2 + kv = FC y) exp( SE 2 ds 
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Le résultat est donc établi, quitte à poser 
ax + By \ 
G(XYN=F(XN exp |. 9 
Généralisation. | 
Soit une équation aux dérivées partielles linéaire à coefficients 


constants à n variables, qui s'écrit donc sous la forme suivante 
(avec, pour tout couple (j, j), a; = a;;): 


+ E ou 
2 2° Frs D de ee (3) 
i= 1 j = 


ne par Q la forme quadratique 


n n 
> Dax Xpr 


i=1 j=1 


DER A) © 


on dira alors que cette équation (3) est: 


— hyperbolique, si sa signature est de la forme (p, n-p) avec 
0O<p<n; 

— parabolique, si sa signature est de la forme (p, q) avec 
p+q<n; 

— elliptique, enfin, si sa signature est (n, 0) ou (0, n). 

(On rappelle que si (p, q) est la signature d’une forme quadra- 
tique, p et q sont respectivement les nombres des carrés précédés 
de signes + et —- dans une décomposition en carrés de formes 
linéaires indépendantes). 


On aborde maintenant quelques exemples illustrant ces trois 
types d’'équation, l'équation des ondes, l'équation de la chaleur et 
l'équation de Laplace. 


Les paragraphes 2, 3 et 4 sont effectivement consacrés à l'étude 
d'équations importantes sur le plan physique, l’ouvert U étant à 
chaque fois géométriquement simple afin de conduire les calculs à 
leur terme. 


squation 
PA: 


lion des ondes 


2.1 Équation des cordes vibrantes 


Il s'agit de l’une des premières équations aux dérivées partielles 
mises en évidence. 


Elle fut étudiée dès la première moitié du XVII® siècle par 
d'Alembert : 


ae" 3% 


où v désigne la vitesse de propagation de l'onde dans la corde et 
u (x, t) l'ordonnée du point d’abscisse x de la corde à l'instant t 
(cette ordonnée étant mesurée par rapport à la position d'équilibre 
supposée d’'ordonnée nulle). 


2.1.1 Cas d'une corde infinie 


# On suppose la corde vibrante infinie, et on assimile la position 
d'équilibre de celle-ci à la droite réelle R. On se propose d'étudier 


l'équation avec les conditions initiales suivantes, supposées réali- 
sées pour tout nombre réel x: 


u(x, 0)=f(x) 
| ou 
et 3e C6: 0) = 0. 


Ces conditions signifient que la corde a été lâchée sans vitesse 
initiale à partir d’une position définie par la donnée de la fonction f, 
que l'on suppose de classe C? sur R (ou même de classe C ! et de 
classe C? par morceaux). 


En s'inspirant des changements de bases qui ont été réalisés au 
paragraphe 1, on considère les nouvelles variables : 


X=x+vtet T=x-vt 
et l’on pose : 
VIX T)=vix+vt, x-vt)=u(xt). 


Æ Par un simple calcul de dérivées partielles, on constate que u 
satisfait l'équation des cordes vibrantes si, et seulement si, v vérifie 
l'équation : 


Il en résulte que ok ne dépend pas de X, et donc que : 


oT 


av 
ST (X, T) = a(T) 


où a est une fonction de classe C1. 


Quitte à noter À une primitive de a, on obtient alors par une 
nouvelle intégration : 


v(X T)=A(T)+B(X) 
où À et B sont donc deux fonctions de classe C? de R dans R. 
On obtient, par conséquent : 
u(x, t)=A(x-vt)+B{(x+vt). 


Réciproquement, il est immédiat de constater qu'une telle fonc- 
tion est bien solution de l'équation des cordes vibrantes. 


Tenant maintenant compte des deux conditions initiales, il vient : 
A(x)+B(x)= f(x) 
et A'{x)-B'{(x)=0 
Il en résulte que 


LA 


09 = BG = 09, 


ce qui montre enfin que : 
2A(x) = f(xX)+21 et 2B(x) = f(x) -À 


et 2B(x) = f(x)-2 
où À désigne une constante réelle. 


Ainsi, le problème posé admet une solution et une seule, donnée 
par: 


u(x, t) = ue vt})+ {x+vt)). 


Cette solution possède la remarquable interprétation physi- 
que suivante. 

e Considérons uniquement le premier terme f(x- vt)}/2. Celui-ci 
donne l'amplitude du mouvement d'une corde à l'instant t et à l’abs- 
cisse x, et l’on remarque que c’est la même amplitude que celle figu- 
rant à l'instant initial 0 à l'abscisse x — vt. 
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ns a ———— ES 
D / # X 
La corde à l'instant initial t = 0 (l’ordonnée à l'abscisse x est f(x)) 
(a) corde à l'instant initial 
Pa Pa 
ré \/ \N 
= — — — = eee 
' pe N # X 
# A N 4 
v 
Cheminement de l'onde vers la gauche aux instants f,, puis t, 
(b) propagation de la première onde 
M PAN 
NN 
PR de “nn Go prima" jé 
b N, } à #7 ù | X X 
VU 4  F V7 \ PF d 
Cheminement de l'onde vers la droite aux instants f,, puis t, 
(©) propagation de la deuxième onde 
En ET  — a 
>» ue, dd SNS > ae. Ff LS à % X 


Superposition des deux ondes aux instants t,, puis t, 
(l'ordonnée à l’abscisse x est (f(x vt) + f(x + vt)) 


() superposition des deux ondes 


Figure 1 - Problème de la corde vibrante infinie 


Ainsi, ce terme indique le cheminement d'une onde qui se 
déplace vers la droite à la vitesse v, puisque, entre les instants 0 et t, 
elle est passée de l’abscisse x - vt à l'abscisse x. 


e De même, le second terme f{x+ vt)}/2 indique le cheminement 
d’une onde qui se déplace vers la gauche à la même vitesse v. 


Le problème de la corde vibrante infinie est représenté figure 1. 


2.1.2 Cas d’une corde finie 


2.1.2.1 Présentation 


On suppose la corde vibrante finie de longueur L, et on assimile la 
position d'équilibre de celle-ci au segment [0, L]. On se propose ici 
d'étudier l'équation avec : 

— les conditions initiales suivantes, réalisées pour tout nombre 
réel x: 


u(x, 0) = f(x) 
au 
et ÉT: (x,0) = 0 


qui signifient que la corde a été lâchée sans vitesse initiale à partir 
d'une position définie par la donnée de la fonction 7, que l'on 
suppose de classe C? sur [0, L] (ou même de classe Cl et de classe 
C< par morceaux) ; 
— les conditions aux limites suivantes, réalisées pour tout nom- 
bre réel positif t: Ù 
u(0,t)=u(L,t)=0 
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qui signifient que la corde est fixée à ses deux extrémités. 
Pour t = 0, on a donc: 


f (0) = f(L) = 0. 


On étudie les problèmes d'existence et d'unicité de la solution u 
d’une telle e.d.p. 


2.1.2.2 Unicité d'une éventuelle solution 
par considération de l'énergie 


Supposons que l’on dispose de deux solutions u; et u, de l'équa- 
tion (qui peut ici être homogène ou avec second membre, ce qui ne 
change rien à la démonstration), vérifiant les mêmes conditions 
initiales et les mêmes conditions aux limites. 


Posons u = u, — u, et introduisons la fonction d'énergie suivante : 


E(t) = ; LE CS (X, D) + (92 (X, D) }ax. 


Les hypothèses faites autorisent la dérivation de E sous le signe 
intégral, et on a, puisque u = u, — u, est aussi solution de l'équation 
des cordes vibrantes : 


; L/1 du d2u au d2u 

ED LE #06 002 0e 0 +8 06 D PE 06 }x 
Lau o2u ou d2u 

: LE CDS a O6 D +5 06 Dax 0 D Ja. 


L'expression sous le signe intégral est une dérivée et on a 
ou au 
ETS t) — Se t) D 0 


par dérivation des relations u (0,t)=u{(L, t) = 0, ce qui donne: 


, ou ou : 
E’(t) = É 06 DS C6 o| be 

Par conséquent, la fonction t— E(t) est constante et, en fait, nulle 
puisque E (0) = 0. Cela résulte de l'expression de E (0) et du fait que, 
u, et u) vérifiant les mêmes conditions initiales, leur différence u 

vérifie : 
ou 
"at 
— d'autre part, u(x, 0) = 0 


— d’une part (x, 0) = 0 ; 
TU au 
donc, par dérivation, sn (x, 0) = 0. 


Puisque E (ft) = 0, on a donc, pour tz0 et0<x<« L : 


La fonction u est par conséquent constante sur [0, L] x [0, +], et 
en fait nulle puisque u (x, 0) = 0. Ainsi, on a bien u, = u, et l'unicité 
annoncée. 


2.1.2.3 Existence d'une solution par la méthode 
de séparation des variables 


Indiquons tout d’abord l'idée de la méthode. 


E On commence par rechercher des solutions multiplicatives non 
nulles de la forme 


(x, t)— u(x) v(t) 


qui vérifient les conditions aux limites vues paragraphe 2.1.2.1. 
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lci, de telles solutions vérifient donc : 


utx)v'"{(t)=vlu" (x) vit). . 


Puisque l’on recherche des solutions non nulles, il existe a priori 
des nombres réels x, et t pour lesquels : 


u(x)4#0 et v(t)#0. 
On en déduit, quitte à fixer x = xp, puis t=t0, l'existence de 
constantes 1 et u telles que l'on ait pour t= 0 et 0<x<L : 
u"(x)=Aut(x) ; 
v"(t)=uvit). 
En reportant réciproquement dans l'équation, on voit que, en fait : 
u=v?2 


et pour tenir compte des conditons aux limites, il faut enfin bien sûr 
que : 
u(0)=u(L)=0. 
(On est ainsi amené à résoudre un problème, ici très simple, de 
Sturm-Liouville (8 5); on verra que cela est général dans cette 
méthode de séparation des variables). 


e Sile nombre À est nul, alors : 
u(x)=Ax+8B 


et puisque u (0) = u(L) = 0, on a A = B=0 et u est la solution nulle, 
ce qui est exclu. 
e Si le nombre À est strictement positif, on a : 


u(x) = Ach(4/Àx) + Bsh(/Ax) 


et puisque u (0) = u (L) = 0, on a A= B=0 et u est la solution nulle, 
ce qui est exclu. 
e Sile nombre À est strictement négatif, on a: 


u(x) = À cos(4/-Àx) + B sin(/-A1x) 


puisque u (0) = 0, on a À = 0 et puisque u (L) = 0, on a B=0, donc 
encore u = 0, sauf s’il existe un nombre entier naturel non nul ntel 
que 


-n?2 nr? 


À = 12 


r 


auquels cas on obtient alors les solutions suivantes : 


u(x) = À, sin(? x) ; 


Ë 


v(t) = B, cos(? Fe )+ Csin(? re ). 


& L'équation étant linéaire, les combinaisons linéaires des solutions 
précédentes sont encore solutions de l'équation. L'idée est de ne 
pas se borner nécessairement à des «combinaisons linéaires 
finies » pour obtenir une solution. 


Posons donc, de façon purement formelle (on peut faire À, = 1): 
._(n ax n Trvt . (nTrvt 
u(xX À = Z sin EC nl )+ Cusin( nl ) : 


Les conditions de Cauchy portant sur u (x, 0) et x 0) seront 


formellement vérifiées en choisissant les coefficients C, nuls et les 
coefficients B, tels que : 


f(x) = ÿB, sin( 2 ). 


Un tel développement est celui d’une fonction impaire et 21L- 
périodique sur R, que l’on obtient en prolongeant la fonction f par 
imparité sur [- L, L], puis par 2L-périodicité. La fonction f ainsi 
prolongée est clairement de classe C sur R. Elle est a priori de 
classe C? sur R privé de l’ensemble LZ des multiples de L, sauf si 
elle vérifie la conditon supplémentaire 


f'A\0=T" =D, 


auquel cas elle est.alors de classe C?surR. 


On en déduit qu’elle est développables en série de Fourier et que 
sa série de Fourier converge normalement vers f. 


Si l'on pose pour t=0 et 0<x<L : 


+ 0° 


u(x, t ) = +. B, sin(? FX jcos[2 7) . 


n=1 


L 
B, = 1 [ f(X) sin(? FX )dx 


on voit que la série définissant u est normalement convergente sur 
[0, L]X[0,+c], donc continue, et l’on vérifie aisément, à l’aide des 
formules de trigonométrie et du développement en série de Fourier 
de f, que l’on a: 


LfCX + VE) + f(x vb] 


u(x, t) = 2 


La fonction u obtenue ci-dessus est donc de classe C? sur 
[0, LIX[0,+c] et vérifie les conditions aux limites et les conditions 
initiales requises. 

Elle est de plus de classe C? sur [0, L]x[0,+c] et solution de 
l'équation des cordes vibrantes si : 

f" (0) =f"(L)= 0. 


Sinon, elle n’est de classe C? et solution de l'équation des cordes 
vibrantes que sur l’ensemble [0, L] x [0,+c] privé des segments de 
droite d'équations 


xtvt= KL 


avec t=0,0<x<LetkezZ. 


| 

1. La formule obtenue pour une corde finie (après avoir pro- 
longé convenablement par imparité et périodicité la fonction f) 
est exactement l’analogue de la formule obtenue dans le cas 
d'une corde infinie (8 2.1.1). 


2. La solution u ainsi obtenue dépend continüment de la don- 
née initiale f. En effet, on remarque facilement que l'on a: 


ul. <1fl. 


Supposons que la donnée initiale f ne soit connue que par 
une approximation f. 

Notons alors u et u les solutions associées à ces données ini- 
tiales exactes et approchées fet fet 


Au =u-u et Af= F-T 


les erreurs. On a, par différence : 
AU(X, t) = LAF (+ vr}+ AF (e= vo] 


et donc, de même : 


IAu |. = sup [u(x)-u(x)] <IAfIS = sup [f C9 - C9). 
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INTRODUCTION AUX ÉQUATIONS AUX DÉRIVÉES PARTIELLES LINÉAIRES 


2.2 Généralisation : 
l'équation des ondes en dimension n 


Il s’agit de l'équation aux dérivées partielles : 


ot? OX? oxe ‘ 
Ainsi, pour n= 2, on obtient l'équation des membranes vibrantes : 


o2u 2 ( ue) 
— =V + 


me "rap 


où v désigne la vitesse de propagation de l’onde dans la membrane 
et u(x, y, t) la cote du point de coordonnées (x, y) d’une membrane 
vibrante à l'instant t (cette cote étant mesurée par rapport à la posi- 
tion d'équilibre supposée de cote nulle). 


Pour une membrane rectangulaire par exemple, dont la position 
d'équilibre s'’identifie au rectangle R = [0, a]x[0, b] du plan, le 
problème est d'étudier l'équation avec : 


— les conditions initiales suivantes, réalisées pour tout nombre 
réel x : 
u (x, y, 0) = f(x, y) 


au 
et 30% y,0) = 0 


qui signifient que la membrane a été lâchée sans vitesse initiale à 
partir d'une position définie par la donnée de la fonction f, que l'on 
suppose de classe C?sur R; 


— les conditions aux limites suivantes, réalisées pour tout nom- 
bre réel positif t: 


u(x, 0,t)= u(x, b,t) = u(0,y,t)= u(a,y,t)= 0 
(O0<x<a,0<y<b). 

qui signifient que la membrane est fixée aux bords du rectangle FR. 

Pour t= 0, on a donc: 

f(x 0) = f(x, b) = f(0, y) = f(a, y) = 0 (O<x<a, 0<y<b) 

L'unicité des solutions de ce problème se traite exactement 

comme au paragraphe 2.1.2.2. supposons que l’on dispose de deux 

solutions u, et u, de |” équation (qui peut ici être homogène ou avec 


second membre, ce qui ne change rien à la démonstration), vérifiant 
les mêmes conditions initiales et les mêmes conditions aux limites. 


Posons u = u, — u, et introduisons la fonction d'énergie suivante : 


E(t) = 3 [LC 2 (Sice y, D) + + (SEcx y, p) + (ice y, D) Jax dy. 


ILest licite de dériver sous le signe intégral et, eñ omettant d'écrire 
(x, y, t) sous le signe intégral, il vient donc : 


es CLS du d2u ou o?u 
E"(® = A À dt 32 9x 3tox" dy ot Je) à 


Puisque u = u, - u, est aussi solution de l'équation des 
membranes vibrantes, on a: 


: (+ d2u\ du d2u du d2u 
ui: dE JC Gt pe) * dx EDR 97 DD) 


t 
| La qua). 2 (Ru 20 ax 
à (5 ot ox ay Kat dy Y- 
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Comme u (0, y, t) = u (a, y, t)= 0, on a, par dérivation : 
ou ou 
ÉTA l'A t) En se ÿ; t) on 0, 


et de même 


ou au 
Fri Ge 0, t) = 31% b, t) æ 0, 


ce qui conduit aussitôt à : 


: b ru au? du au? 
FAT L (5: ar (3 Sy, bu 
Par conséquent, la fonction t— E(t) est constante et, en fait, nulle 


puisque E (0) = 0. Cela résulte de l'expression de E (0) et du fait que 
u, et u, vérifiant les mêmes conditions initiales, on a : 


— d'une part, _ (X, y, 0) = 0 ; 
— d'autre part, u (x, y, 0) = 0 
donc par dérivation : 


au ou 
3x 06 Y, 0) = 0 et y V2 0) = 0. 


Puisque E (t) = 0, on a donc, pour t= 0 et0<x<a,0<y< pb : 
ou ou ou 
spa 2 t) = ay 2 t) _ ÉTIOS y; t) = 0. 


La fonction u est par conséquent constante sur 
[0, a]x[0, b] X[0, +cel, 
et elle est nulle puisque u (x, y, 0) = 0. Ainsi, on a bien u; = u, et 
l’unicité annoncée. 


De même, la méthode de séparation des variables permet 
d'obtenir la solution u du problème lorsque la fonction fest suppo- 
sée suffisamment régulière. La résolution du problème de Sturm- 
Liouville conduit à nouveau aux fonctions trigonométriques, alors 
que, avec une membrane circulaire, elle conduirait aux fonctions de 
Bessel. 


3. Une équation 
l'équation de 


3.1 Équation de la chaleur en dimension 1 


Il s'agit de l'équation aux dérivées partielles suivante : 


Les constantes c et ydésignent la capacité et la conductivité ther- 
miques d’une barre et u(x, t) indique la température au point. 
d'abscisse x de la barre à l'instant f. 


3.1.1 Cas d'une barre infinie 


On suppose la barre infinie, et on assimile celle-ci à la droite réelle 
R. On se propose d'étudier l'équation avec pour condition initiale, 
réalisée pour tout nombre réel x: 


u (x, 0) = f(x). 
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Celle-ci indique la température des différents points de la barre à 
l'instant initial. On supposera cette fonction fcontinue et bornée sur 
R et, pour simplifier les notations, on posera 


Hs 


C 


Pour trouver la forme d'une solution, nous allons suivre ici une 
méthode différente et plus complexe que celle mise en œuvre dans 
le cas de l'équation hyperbolique correspondante (équation des cor- 
des vibrantes (8 2.1)), reposant sur la transformation de Fourier, 
à propos de laquelle nous faisons quelques brefs rappels. 


e La transformée de Fourier Ff d'une fonction sommable f est 
définie par : 


FF (x) = [' © f (Yexp(-2irxy)dy. 


e La transformée de Fourier conjuguée de la fonction Ffsupposée 
sommable donne : 


f (x) = (à © Ff (Yexp(2irxy)dy. 


e Par un calcul très classique, on établit que la transformée de 
Fourier de la fonction x— exp(-ax2) est la fonction 


T x? x? 
x E exp à ) (a>0). 


# De plus, nous allons faire des hypothèses plus fortes en suppo- 
sant, d’une part, f sommable sur R, d'autre part, en cherchant des 


solutions u sommables sur R ainsi que leurs dérivées et telles que 


du soit dominée sur R par une fonction sommable U indépendante 


ot 
de t. 


# En appliquant la transformation de Fourier à l'équation de la cha- 
leur par rapport à la variable x, on obtient alors : 


co co 2 
"au (y, texp(—-2irxy)dy = 2 us = (y, ©) exp(-2irxy)dy. 
SD | 4% 


La condition de domination E (y, texp(—-2irnxy) = U(y) par la 


fonction sommable Ü permet de justifier, au premier membre, la 
permuation des symboles d'intégration et de dérivation par rapport 
à la variable t, tandis que la sommabilité de u et de ses deux 
premières dérivées en x permet d'appliquer deux fois, au second 
membre, la formule F(f’)(x%) = 2irxFf(x), permettant de simpli- 
fier la formule précédente qui devient : 


à f+e 
Si u (y, exp(-2irxy)dy 


ne u (y, texp(-2irxy)dy. 


= an2o x Ï 

Notant v la transformée de Fourier de u par rapport à sa première 
variable, on a donc: 
dv 

ot 


Cette équation différentielle se résout aisément et donne : 


(x, t)+ 4n2@2 x2 V(X, t) = 0. 


v(x t) = Ff(x) exp(- 4n2o2 x2t) 


puisque, pour t= 0, l'expression de droite v (x, 0) n’est autre que la 
transformée de Fourier de la fonction x— u(x, 0) = f(x). 
On note alors que : 


— la fonction x v(x, ft) est sommable (pour t > 0) comme pro- 


duit de la fonction bornée Ff (puisqu'une transformée de Fourier 
tend vers 0 en ++, donc est bornée) par la fonction sommable 
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x exp(- 4x2w2x2t) ; par transformée de Fourier inverse (ou 
conjuguée), on a donc l'expression de u (x, t} à partir de celle de 
vix t)sit>0; | 

— la fonction x— exp(- 4r2w2?x?t) est, d'après les rappels pré- 
cédents, la transformée de Fourier de la fonction 


2 2 
cé exp(-x*/4@ Ù. 
Nana? t 


Dans ces conditions, la transformée de Fourier (directe ou conju- 
guée) du produit est le produit de convolution des transformées de 
Fourier (directes ou conjuguées), d’où : 


1 Li _ @- y} 
20.fri M ne SE Jar 


On notera que l'application de la transformation de Fourier peut 
s'effectuer avec bien d'autres équations aux dérivées partielles, 
mais elle ne trouve son véritable champ d'application que dans le 
cadre des distributions. 


X 


u(x, t) = 


& La dernière formule ci-dessus définit en fait u (x, t) pour tout cou- 
ple (x, t) de nombres réels tels que t > 0 dès que fest supposée seu- 
lement continue et bornée sur R. 


Il s'agit maintenant de savoir, réciproquement, si cette fonction u 
est de classe C? et si elle est solution de l'équation de la chaleur et 
si, enfin, u (x, t) tend bien vers f(x) lorsque (x, t) tend vers (x, 0). 

e Le théorème de dérivation des intégrales dépendant d’un para- 
mètre montre que u est de classe C? sur Rx] 0,+  [ et les déri- 
vées de l'intégrale s’obtiennent par dérivation sous le signe intégral, 
ce qui permet de vérifier par un simple calcul que u est solution de 
l'équation de la chaleur. 


e La fonction u vérifie la condition de Cauchy désirée, car u (x, t) 
tend vers f(x) lorsque (x, ft) tend vers (xo, 0). 


En posant en effet x- y = 2 zJt ,ona: 
u(X, t) = E ['Trcc-20 zNt )exp(-22 )dz. 
x J-= 


Sous cette dernière forme, le théorème de continuité des inté- 
grales dépendant d’un paramètre prouve immédiatement le résultat 
voulu. 


Unicité des solutions 


M Si u, et u, sont deux solutions de notre problème, leur diffé- 
rence u = u, — u, est solution de l'équation de la chaleur avec 
u (x, 0) = 0. 


& Si u, et u, sont sommables sur R ainsi que leurs dérivées et 


tell au; ë ou, 
elles que — et —< 

M Se Ÿ 
sommables indépendantes de t, les calculs menés précédem- 
ment montrent que la transformée de Fourier de u est nulle, 
donc que uest nulle, ce qui établit - avec ces hypothèses supplé- 


mentaires — l’unicité de la solution du problème posé. 


Dans le cas général, on ne peut a priori conclure à l’unicité des 
solutions. Ainsi, la fonction définie par : 


soient dominées sur R par des fonctions 


#5 x? 
Mée LT- T7 


pour t> 0 est bien solution de l'équation de la chaleur, et on peut 
vérifier que u(x, t) se prolonge par continuité pour {= 0 par 
u (x, 0) = 0. Mais la fonction nulle est aussi solution de ce même 
problème de Cauchy, pour lequel il n'y a donc pas a priori unicité 
de la solution sans hypothèses supplémentaires (portant ici sur 
la croissance des solutions considérées). 
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3.1.2 Cas d’une barre finie 


On suppose la barre finie de longueur L, et on assimile celle-ci au 
segment [0, L ]. On se propose ici d'étudier l'équation avec : 


— les conditions aux limites : 
u(0,t)=ut(L,t)=0 
pour tout nombre réel positif {, qui indiquent que les extrémités de 
la barre sont maintenues à la température 0 ; 
— la condition initiale : 
u (x, 0) = f(x) 
pour 0< x< L, qui indique la température des différents points de 
la barre à l'instant initial ; 
— avec aussi, bien sûr : 
f(0)=f(L)= 
On supposera ici cette fonction f continue et de classe C! par 


morceaux sur [0, L] et, pour simplifier les notations, on posera 
encore 


st. 
c 


3.1.2.1 Unicité d'une éventuelle solution 


Deux méthodes peuvent ici être mises en œuvre, l’une analogue à 
celle développée dans l'étude des équations hyperboliques (8 2.1.2.2), 
l’autre analogue à celle qui sera développée pour les équations 
elliptiques (8 4). 


On suppose donc que l'on a deux solutions u, et u, de |” équation 
(qui peut ici être homogène ou avec second membre, ce qui ne 
change rien à la démonstration), vérifiant la même condition initiale 
et les mêmes conditions aux limites. 


On pose alors u = u, — u, et l'on a donc: 
u(x,0)=0etu(0,t)=ui(L,t)=0 


M Première méthode 
On introduit la fonction suivante, obtenue en multipliant l'équa- 
tion de la chaleur par u et en l’intégrant par rapport à x sur [0, L]: 


Le ou (X, t) u(x, t)dx = œ? É 2 g (x, t) u(x, dx. 
0 0 


Le premier membre apparaît comme une dérivée (on peut ici 
permuter les symboles de dérivation et d'intégration) et l'on peut 
intégrer par parties le second membre en tenant compte du fait que 
u(0,t)=ui{L,t)=0,ce qui donne: 


L 
12 u? Ge des er [| (2 660 ÿ ax x<0. 
0 0 0 


Cela prouve que la fonction e [| u? (x, t)dx est décroissante. 


Comme elle est nulle en 0 et visiblement à valeurs positives, elle est 
donc nulle, ce qui implique la nullité de u et l'égalité u = u2.. 


& Seconde méthode 
On va utiliser ici le principe du maximum. 
Considérons un nombre réel strictement positif £ et la fonction : 
u, (xt) = u(x )+ex?, 
qui vérifie : 
2 


Le 
ax? 


ue 2 d Us L ou 
ot (x t)- © (x, t) à TAG t)- © 


5e (x. t)-2eo? = _2e0?. 
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Sur le compact [0, L]x[0, T ] où T > 0, elle admet un maximum, 
atteint en un point (x, t)) dont nous allons montrer qu'il est tel que : 


x) = Dou LavecO<t,<T 
ou tel que : 
= Oavec0<x,<L. 
Si tel n'est pas le cas, on a en effet : 


0O<x,<L et O<t<T 


et: 


— la fonction x—u,(x t,) atteint son maximum en x sur 
l'intervalle ouvert 10, LI de sorte que l'on a les deux relations : 


3x Po fo) = 0 ; 


— la fonction t—u,(x,t) atteint son maximum en t sur 
l'intervalle ]0, T[ de sorte que, en considérant sa dérivée en t 
comme sa dérivée à gauche, on a par définition : 


On a alors, en combinant ces deux inégalités : 


2 


9° Ue 7. 


Cette contradiction prouve le résultat annoncé, à savoir que U: 
atteint son maximum nécessairement sur l'un des trois bords infé- 
rieurs du rectangle [0,L]x[0, T]. Comme on sait que 
u, (X ©) = u(x t)+ex? et que u est nulle sur les trois bords infé- 
rieurs de ce rectangle, ce maximum de u, est donc inférieur ou égal 
à eL*. 


Il en résulte que le maximum de u sur ce même rectangle, qui est 
inférieur ou égal à celui de u, est lui-même inférieur à EL?, Comme € 
est arbitraire, ce maximum de u est donc négatif et u est donc à 
valeurs négatives sur le rectangle [0, L]x[0, T7 ]. 


Par le même raisonnement, quitte à remplacer u par — u, on voit 
que - u est aussi à valeurs négatives, et donc u à valeurs positives, 
sur ce même rectangle. 


En conclusion, pour tout T=0, u est nulle sur [0, L]x[0,T], 
donc sur [0, L]x[0, + c ] et l’on a bien l'égalité u4 = u. 


3.1.2.2 Existence d'une solution 
par la méthode de séparation des variables 
L'idée de la méthode est la même que pour les équations hyper- 
boliques (8 2.1.2.3). 


On commence par rechercher des solutions multiplicatives non 
nulles de la forme 


(x D) u(x) v(E) 


vérifiant les conditions aux limites vues au début du para- 
graphe 3.1.2. 


lci, de telles solutions vérifient donc : 


u (x) v'(t)=œ2u"(x) vit). 


Puisque l’on recherche des solutions non nulles, il existe a priori 
des nombres réels x et f, pour lesquels : 


u(X%)40 et v(tH)#0. 
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On en déduit, quitte à fixer x= xp, puis t= tp, l'existence de 
constantes 1 et u telles que l’on aït pour t= 0 et 0< x< L : 


u"{(x)=Au(x) ; 
v'(t)=uvi(t). 
En reportant réciproquement dans l'équation, on voit que, en fait : 
u = ©? À 


et pour tenir compte des conditions aux limites, il faut enfin bien sûr 
que : 


u(0)=u(L)=0. 


(On est encore amené à résoudre un problème, très simple, de 
Sturm-Liouville). Cette dernière condition ne peut être évidemment 
réalisée (ainsi qu'on l’a déjà détaillé dans le paragraphe 2.1.2.3) que 
si, et seulement si : 


_— nr? 


À = 12 


où n désigne un nombre entier quelconque supérieur ou égal à 1, et 
l’on obtient alors les solutions suivantes : 


u(x) = À, sin( TX ) 1 


2 hn2 r2 
— On Tr< t 
v(t) = B, exp(— ). 
 L'équation étant linéaire, les combinaisons linéaires des solutions 
précédentes sont encore solutions de l'équation. L'idée est à nou- 
veau de ne pas se borner à des « combinaisons linéaires finies » 
pour obtenir une solution. 


Posons donc, de façon purement formelle (on peut faire A,;= 1): 


2 nn 5 
u(x t) = YB, sin[ EX Jexp|- + ). 


La condition de Cauchy portant sur u{(x, 0) sera formellement 
vérifiée en choisissant les coefficients B,, tels que : 


fx) = YB, sin[ 2x ). 


Un tel développment est celui d'une fonction impaire et 2L-pério- 
dique sur R , que l’on obtient en prolongeant la fonction fpar impa- 
rité sur [- L, L ], puis par 2L-périodicité. La fonction fainsi prolongée 
est clairement continue, de classe C! par morceaux sur R. On en 
déduit qu'elle est développable en série de Fourier et que sa série de 
Fourier converge normalement vers f. 


Si l’on pose donc : 


< . (NTX @? n?r2t 
u(xt)= Ÿ B, sin 1 JP — Je 
n=1 


2 . NY 
B,= 1 [ro Sin dy, 


on voit que la série définissant u est normalement convergente sur 
[0, L] X[0, + c ], donc continue, et l’on vérifie aisément, à l’aide des 
théorèmes sur les séries de fonctions, qu'une telle fonction est de 
classe C? sur [0, L] x [0, +  ], les dérivations s’effectuant terme à 
terme, ce qui permet de vérifier que la fonction obtenue, qui vérifie 
les conditions initiale et aux limites, est solution de l'équation de la 
chaleur. 


On notera que la convergence normale de la série, figurant ci- 
dessus sous le signe intégral, autorise la permutation des symboles 
de sommation et d'intégration, de sorte qu’en utilisant les formules 
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de trigonométrie usuelle, on obtient (la fonction f étant toujours 
supposée prolongée sur R par imparité et 2L-périodicité) les 
égalités suivantes : 


: + L Fit . 202 
u(Xx, t) = sl Î Ré cos UN y) exp(- MER * Ja 
S n=1 


L [2 


PRE Ne ) 5 ex inm(x= y) ex (- D Va 
= 2L nr ee [e) L p L2 


(Pour n=0, l'intégrale de la fonction impaire fsur [- L, L] est nulle. 


Cela s'écrit également avec la formule sommatoire de Poisson 
(cf. encadré) : 


1 *L a (x-y-2nL)2 
RD = f ex (CRE A 
ET lé D 2 F Lot à 


lci encore, la convergence normale de la série figurant sous le 
signe intégral justifie la permutation des symboles d'intégration et 
de sommation. 


En effectuant alors, dans chaque intégrale, le changement de 
variables 
z=y-2nL 


puis en prenant en compte la 2L-périodicité de f (dont on rappelle 
qu'elle a été convenablement prolongée sur R auparavant), on 
obtient enfin la formule suivante : 


1 +. _ (x-2z} 
20 nt IRC exp| 4 t je 


On remarquera qu'elle est très exactement l’analogue de celle 
obtenue dans le cas d’une barre infinie (8 3.1.1). 


u(x ?) = 


1. La formule sommatoire de Poisson que l’on vient d'utiliser 
s'obtient facilement par développement en série de Fourier de 
la fonction indéfiniment dérivable et 2L-périodique définie avec 
a>0par: 


+ 00 


g(u)= Y exp(-a(u-2nLy?). 


n=-—0 


Le calcul de son coefficient de Fourier c, donne en effet : 


1 fr£ ikru 
Ck = 5L _ exp ee) g(u)du 


1 te 2 ikru 
“x |. + exp(-a(u-2nL)?2) exp{- FU Jau. 


NN =—-0c 


Par convergence uniforme, il est possible de permuter les 
symboles d'intégration et de sommation, ce qui donne alors : 


4 << fe ane ikru \, 
= 3] + _, xP(atu- nL}?) exp|- u 


N=—0 


En effectuant le changement de variable v=u-2nlL, on 
obtient (en tenant compte de la périodicité de l’exponentielle 
complexe) : 


+ oo 
1 + (2n+1)L à ikrv 
C = — exp(-av<) ex (av 
nr À nee de So: 


EL ds 2 IKnv 
= 35 [ expt-ave) exp(- L Ja. 
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Quitte à intégrer la forme différentielle (clairement exacte) 
exp(- az?)dz sur le bord du rectangle dont les quatre sommets 
sont (- À, 0), (A, 0), (A, ikx/2aL) et (- À, ikn/2aL), puis à faire ten- 
dre À vers + (les deux intégrales sur les bords verticaux du 
rectangle ont alors pour limite 0), on obtient : 


sis ikn \? nu de 2 - 
Le exp{-a[v+ nn) Jav = lé exp(-av< )dv = si 
Il en résulte que la valeur du coefficient de Fourier cest égale 
à : 
CL = 1 ex (- L E) "Tex (at ) Jav 
CAR: dos FT ES D Etes 2aL 


“x LÉ 
dE P 4aL?2 /JNa° 


Puisque la fonction à est développable en série de Fourier 


‘ d’après le théorème de Dirichlet, on obtient la version suivante 


de la formule sommatoire de Poisson : 


+ co 


+ © 
>, exp(-a(u-2nLy?) = RTE ÿ exp 
k=-—0 


Nn=—-c 


(- k2 n2 ) ae (es 
4aL? L 
et la forme utilisée dans le calcul mené plus haut en résulte en 

faisantu=x-yeta= 1/40? t. 

2. Contrairement à la situation qui prévalait avec l'équation 
des ondes (de type hyperbolique), on notera qu'avec l'équation 
de la chaleur (de type parabolique), il est hors de question de 
vouloir « remonter le temps » en attribuant à t des valeurs néga- 
tives, car les résultats obtenus perdent alors tout sens. 


3. La solution u ainsi obtenue dépend continûüment de la don- 
née initiale f. En effet, on remarque facilement que l'on a: 


Iul. <1Fl. 


Comme on l’a déjà indiqué dans le cadre de l'équation des 
cordes vibrantes, supposons que la donnée initiale f ne soit 
connue que par une approximation f. Notons alors u et u les 
solutions associées à ces données initiales exactes et appro- 
chées fet fet 


Au=u-uet Af=f-f 


les erreurs. On a par différence : 


1 fee 
Af _ 
 . F aexp( 


Et donc de même 


Au. 


Au (x, f) = CD Jar. 


4? t 


= sup [u(x)-u(x) < 


0<x<L 


IAf[.. = sup |fC9- FC9)|. 


0O<x<L 


3.2 Généralisation : équation 


de la chaleur en dimension n 


ll s’agit de l’équation aux dérivées partielles suivante : 


ou _ Le re) 
0x? da À. 
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| 


) 


# Ainsi, pour une plaque rectangulaire que l’on identifie au rec- 
tangle À = [0,a]x[0, b], le problème est par exemple d'étudier 
l'équation : 


au 7 2 d2u ) 
— = + | — + — 
oxX2 dy? 
où u(x, y, t) désigne la température au point des coordonnées 
(x, y) de la plaque à l'instant t avec : 
— les conditions aux limites suivantes, réalisées pour tout nom- 
bre réel positif t: 


u(x, 0, t) = u(x, b,t) = u(0, y, t) = u(a, y,t)= 0 


(0<x<a0<y< D) 


qui signifient que les bords de la plaque sont maintenues à la tem- 
pérature 0. 

— la condition initiale suivante, réalisée pour tout couple (x, y) 
de nombres réels tels que 0<x<aet0<y<b: 


u (x, y, 0) = f(x, y) 


qui indique la température des différents points de la plaque à l'ins- 
tant initial ; 


— avec aussi, bien sûr : 
f(x 0)= f(x, b)=f(0,y)=f(a, y)=0 


On supposera ici cette fonction f de classe C' sur R et, pour 


simplifier les notations, on posera encore @? = T 


# l'unicité des solutions de ce problème se traite exactement 
comme au paragraphe 3.1.2.1. Supposons que l'on dispose de deux 
solutions u et u, de l'équation (qui peut ici être homogène ou avec 
second membre, ce qui ne change rien à la démonstration), vérifiant 
la même condition initiale et les mêmes conditions aux limites. 


On pose u = u; — u, et l'on a donc: 
u(x 0,t)=u(x b,t)=u(0,y,t)=u(a y,t)=0 
et 
u(x y,0) = 0 (0<x<a,0<y<b). 


On peut alors utiliser l’une des deux méthodes déjà indiquées. 
a) On considère la fonction définie pour t= 0 par: 


JL (x, y, t) U(X, y, dxdy 


[Les 


Le premier membre apparaît comme une dérivée (on peut 
permuter les symboles de dérivation et d'intégration) et l'on peut 
intégrer par parties le second membre en tenant compte du fait que 
u s'annule sur le bord du rectangle ce qui donne: 


2 
55 ! Î[. u? (x, y, t) dxdy 


= — IA 3x tx o) i ) +(5cx o) Jaxdy < 0 ; 


Cela prouve que la fonction [| u2 (x, y, t) dxdy est décrois- 
R 


= (% y, t) +5 — — (X y, Jucx y, t)dxdy. 


sante. Comme elle est nulle en 0 et à valeurs positives, elle est donc 
nulle, ce qui implique la nullité de u et l'égalité u; = u. 
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b) On utilise le principe du maximum en considérant un nombre 
réel strictement positif e et la fonction 


u.(x y,t)=u(x y,t)+e 0 + y?) 
qui vérifie : 
d2u, 


ou, 
= _@? —— 
er (x, y, Ô-@ PE 


2 d7u 2 
(x Lt})-&@ 5,2 (x y, Ü = — 4e@ 
0Y 
Sur le compact [0,4]x[0, b]x[0, T] où T>0, elle admet un 
maximum, atteint en un point (x, Yo, to). On se propose d'établir 
qu'il est tel que : 


x, = Oou aou y, = 0 ou bavecO<t,<T 
ou tel que : 
ty = Dave 0<x, <aet0<y <= b. 


Sinon, on a en effet 0 < x) <aet0<yy<bet 0<t,<T et: 


— la fonction x u, (x Yo tp) atteint son maximum en x, sur 
l'intervalle ouvert ]0, af, la fonction y — u, (x, y, to) atteint son 
maximum en y, Sur l'intervalle ouvert 10, b[, de sorte que l’on a les 
relations : 


du, 
3x C0: Vo: to) = 0; 


d2u 


52 Vo: Vo: to) <0 ; 


au, 
37 V0: Yo: to) = 0 ; 


o2u, 
ET (X0 > Yo» to) <0 ; 
— la fonction tu, (x, Yo t) atteint son maximum en t, sur 
l'intervalle ]0, TI de sorte que, en considérant sa dérivée en t 
comme sa dérivée à gauche, on a par définition : 


du, 
3 Co: Yo’) 7 0. 


On obtient alors en combinant ces deux inégalités : 


du, ; 0° 
ET (Xo » Yo » To) — © 5e (X0 » Yo » to) 


; o2u, s 
Cette contradiction prouve le résultat annoncé, à savoir que u, 
atteint son maximum sur l’un des cinq côtés inférieurs du parallélé- 
pipède [0, a] x[0, b] x [0, T ]. Comme on sait que : 


u(X y.) = u(X y, D+e(x2+y2) 


et que u est nulle sur les cinq côtés inférieurs de ce parallélépipède, 
ce maximum de u, est inférieur ou égal à £ (a? + b?). 


Il en résulte que le maximum de u sur ce même parallélépipède, 
qui est inférieur ou égal à celui de u, est lui-même inférieur à 
e(a2 + b2). Comme e est arbitraire, ce maximum de u est donc 
négatif et u est donc à valeurs négatives sur [0, a]x[0, bIx[0,T ]. 


Par le même raisonnement, quitte à remplacer u par — u, on voit 
que - u est aussi à valeurs négatives, et donc u à valeurs positives, 
sur ce même parallélépipède. 


En conclusion, pour tout T>0, u est nulle sur le parallélépipède 
[0, a] X[0, B] x[0, 7 ] donc sur [0, a]x[0, b] X[0, + +] et l’on a bien 
l'égalité u4 = uo. 
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& De même, la méthode de séparation des variables permet 
d'obtenir la solution u du problème lorsque la fonction fest suppo- 
sée suffisamment régulière. La résolution du problème de Sturm- 
Liouville conduit à nouveau aux fonctions trigonométriques, alors 
qu'avec une plaque circulaire, elle conduirait aux fonctions de Bes- 
sel. 


4.1 Présentation 


Reprenons l'équation de diffusion de la chaleur dans une plaque, 


. que l’on a étudiée paragraphe 3.2. Si l’on parvient à un état d’équi- 


libre, on voit que la température asymptotique d'équilibre u (x, y) 
au point de coordonnées (x, y ) de la plaque vérifie l'équation ellip- 
tique suivante, dite équation de Laplace : 


o2u £ o2u 
ax? ay? 
ou Au= 0. 
Comme on le voit sur cet exemple, cette équation elliptique est 
indépendante du temps et décrit un phénomène stationnaire. 
& De façon générale, en dimension n, l'équation de Laplace s'écrit : 


AU = du, do, 5 du = 
ox? 0x5 ax? 


Dans un premier temps (8 4.2), on étudiera simplement cette 
équation dans un ouvert U (tout en se limitant, pour les démonstra- 
tions, au cas n=2 du plan). Ses solutions sont, par définition, les 
fonctions harmoniques dans l’ouvert Ü et elles jouissent de 
nombreuses propriétés particulièrement intéressantes. 


Dans un second temps (8 4.3) et (8 4.4), on étudiera l'équation de 
Laplace en ajoutant des conditions aux limites à la frontière de 
l’ouvert U (problèmes de Dirichlet ou de Neumann). 


& Exemples de fonctions harmoniques 


Cherchons les fonctions u harmoniques dans R7 sauf à l’origine 
et ne dépendant que de la distance r à l’origine. Celles-ci vérifient 


AÇU(,fx2 + x3 ++ x2)) = 0, 
soit : 


u”’(n+n-lu’(n =0. 


F 
Pour n=2, on obtient : 
u(r)=inr; 
et pour n>3 : 
1 
u(r) = 
(9) (n-2)m-2 


4.2 L'équation de Laplace et les fonctions 
harmoniques dans un ouvert U 


Les solutions de l'équation de Laplace dans U, c'est-à-dire les 
fonctions harmoniques dans U, satisfont plusieurs propriétés 
remarquables, parmi lesquelles les suivantes. 
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Proposition 3. 
On désigne par U un ouvert borné. 


a) Une fonction u est supposée harmonique dans U et se prolon- 
geant par continuité sur l’adhérence Ü de U, c'est-à-dire sur U et sa 
frontière, vérifie le principe du maximum : autrement dit, la fonction 
u atteint son maximum et son minimum sur la frontière de U. 


b) Deux fonctions harmoniques (ou plus généralement vérifiant le 
principe du maximum) dans U et se prolongeant par continuité sur 
l’adhérence Ü de U sont égales dans UÜ si, et seulement si, elles 
sont égales sur la frontière de U. 


Démonstration. © 


Le principe du maximum s'établit comme on l'a fait para- 
graphe 3.1.2.1. On considère un nombre réel strictement positif € et 
la fonction : 


u(x y)=ut(x, y}+e(x?+ y?). 
Puisque u est harmonique, on a: 
AU, = 4€. 


a) Sur le compact Ü constitué de U et de sa frontière, la fonction 
continue u, admet un maximum, atteint en un point Molxo, yo) qui 
appartient nécessairement à la frontière de U. Sinon; (x, yo) appar- 
tient à l’ouvert U et: 


— la fonction x— u, (x yo) atteint son maximum en x sur un 
intervalle ouvert centré en x; 

— la fonction y— u, (x, y) atteint son maximum en x sur un 
intervalle ouvert centré en ÿ,, de sorte que l'on a les relations : 


du, 
ax (X0 ; Yo) = 0 , 


o2u, 
3e (Xo > Yo) <0 ; 


au, 
dy (X0 ; Yo) 0, 


d2u, , sé 
— + (X; <0. 
a Yo 

On en déduit que Au, = 4e<0, ce qui est contradictoire et 
prouve donc que la fonction u, atteint son maximum en un point 
Molxo, Yo) qui appartient comme annoncé à la frontière de U. On a 
donc pour tout point M appartenant à U : 


u(M)< u, (M)< u, (Mo) = u (Mo) +e(xé + 8) u (M5) + ED? 


où D? est la borne supérieure de l’ensemble des carrés des 
distances x? + y? de l’origine O aux points de l’ouvert borné U. 


En faisant tendre € vers 0, on obtient : 
u(M)<u (Mi), 


ce qui établit que u atteint son maximum sur la frontière de. U (on 
procède de même pour le minimum). 


b) Passons au second point. Soient u: et u, deux fonctions harmo- 
niques (ou vérifiant le principe du maximum) sur U qui se prolon- 
gent par continuité sur U. L'égalité de ces deux fonctions sur la 
frontière de U implique que la différence u = u, - u, est harmonique 
dans U et nulle sur la frontière de U : comme u atteint son maximum 
et son minimum sur la frontière de U (principe du maximum) on en 
déduit que u est nulle sur U, et donc que u; = u. 


La réciproque est évidente. () 


Avant d'établir d'autres résultats, on démontre maintenant deux 
formules de Green. 
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Formules de Green : 


=> —> av 
Vu -Vv dxdy+ u Av dxdy = == is . 
D D on 


oD+ 
av au 
u Av dxdy- v Au dxdy = (u Tv Sas. 
D D on on 


oD+ 
Dans ces formules, on désigne : 


— par D un compact connexe du plan dont le bord est supposé 
« suffisamment régulier » ; 
— par dD+ et le bord orienté de D; 


— par n la normale unitaire à 9D+ orientée vers l'extérieur de 
D; 

— par uet v deux fonctions de classe C? sur l'intérieur de D et se 
prolongeant par continuité sur D. 


Démonstration. © 


Ces deux formules résultent directement de la formule de Green- 
Riemann que, selon l'habitude, on écrit sous la forme : 


0Q 9P 
[LGS-S Jaxay = $ Pdx+ Ody. 


0D+ 


La première formule s'obtient en posant ci-dessus : 


av ov 
P=-u— et O=u — 
dy 2 
puis, en permutant les rôles de u et v et en soustrayant les deux 
égalités obtenues, on a la seconde formule. () 


En prenant, dans la première formule de Green, u = 1 et vhar- 
monique, on obtient : 


Cette égalité signifie que le flux du gradient d'une fonction 
harmonique au travers d’un contour (ou d’une surface) fermé 
est nul. 


Proposition 4. 


Une fonction harmonique u dans un ouvert U vérifie le principe de 
la moyenne. Autrement dit, pour tout point Molxo, Yo) de U et pour 
tout nombre r > 0 tel que la boule fermée de centre M, et de rayon r 
soit incluse dans l’ouvert U, le nombre u (x, yo) est la moyenne des 
valeurs prises par u sur la sphère (et aussi sur la boule) de centre M, 
et de rayon r, ce qui, en dimension 2, s'écrit : 


+ 
U(Xo ; Yo) = _ [. u(xQ + rcos 6, y, + rsin6)rd8 


1 . 
“ox (Le RL y) dxdy. 


Réciproquement, on peut établir que toute fonction continue dans 
U qui vérifie le principe de la moyenne est harmonique dans U. 


Démonstration. © 
Dans la seconde formule de Green, convenons de choisir: 


M 
— pour compact D, la couronne £ <M,M<r; avec 0<e<r, 


dont le bord est composé, d’une part du cercle C (M, r) parcourue 
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dans le sens direct, et, d'autre part, du cercle C (M,, £) parcourue 
dans le sens indirect ; 

— pour u, la fonction harmonique de l'énoncé ; . 

— pour v, la fonction (dont on sait qu'elle est harmonique dans 
D) définie par 


v(x, y) = Inr où r= /x2+ y2. 


La seconde formule de Green donne dans ce contexte : 


$ (uv ds - $ (u vos Jds 
on on an on 
0D+ 


CM r )+ 


_ (u op e Jas - 0 
on on 


C(Mo € )+ 


La dérivée de v suivant la normale extérieure à la couronne D en 
un point (x, y) du cercle de centre M, et de rayon r (resp. £) est : 


av dv ] 
— = — = - (resp. 1/E), 
de ar r 


ce qui conduit à l'égalité suivante : 


bL uds-Inr EE ds 
F on 
C(Mo r )+ C(Mo Fr )+ 
A uds +1Ine Es = "0; 
à on 
CM, € )+ C(Mo, € )+ 


Comme le flux du gradient d’une fonction harmonique au travers 
d'un contour fermé est nul (cf. encadré précédent), le premier et le 
troisième des termes figurant dans cette dernière formule sont nuls 
et l’on a enfin: 


_ uds-* uds = 0. 


C(Mo r )+ C(Mo € )+ 


En paramétrant les deux cercles de rayons £et r par l'angle polaire 
6, il vient : 


T7 TH 


Faisant tendre € vers 0 et tenant compte de la continuité de u en 
Mo, le théorème de continuité des intégrales dépendant d’un para- 
mètre donne enfin : 


+T 
Ï u(x, + rcos0, yÿ+rsin0)d6 = 27 u(x, Yo): 


A 


Pour obtenir la seconde égalité, il suffit de changer r en p dans 
l'égalité ci-dessus, de multiplier celle-ci par p et d'intégrer de O à r, 
ce qui donne : 


r Ê+T r 
Ï Î u(xy+pcos8, yo+psin6)pd6 dp = 2 | PU(Xw» Yo)d? : 
0J-7 0 


Passant de coordonnées polaires en coordonnées cartésiennes, 
on a bien: 


 Mdxdy = tr? , : 1 
[Lou 2e y) xdy = Tru (X5 ;, Vo) | 
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1. Le principe de la moyenne apparaît physiquement très 
naturel si l’on se souvient que l'équation de Laplace représente 
des phénomènes stationnaires. 


2. Le principe de la moyenne implique le principe du maxi- 
mum. 


Supposons en effet que la fonction u, harmonique dans un 
ouvert connexe borné U et se prolongeant par continuité sur U, 
atteigne son maximum y dans U. 


Si M, est un point appartenant à U où u atteint son maximum, 
on a: 


u(M)<u(M;) 


sur un Voisinage de ce point M. Choisissant un nombre réel 
strictement positif r assez petit, on a alors : 


1 
— u(x, y)dxdy = u (X, : 
_— JL er 4 PdXdy € U (x: Vo) 

Mais cette inégalité est en fait une égalité d’après le principe 
de la moyenne, et cela implique que, dans cette boule de centre 
Mo et de rayon r,ona: 


u(M) = u(Mi). 


Ainsi, l'ensemble des points M appartenant à U tels que 
u (M) = y est à la fois ouvert (d'après ce qui précède) et fermé 
(d’après la continuité de u). 


Puisque U est connexe, l’ensemble de ces points M tels que 
u (M) = y et égal à U et en fait à ÜU par continuité de u. Ainsi, 
dans ce cas, u est constante. Il en résulte que si u n'est pas 
constante, elle atteint bien son maximum (et son minimum) sur 


+T + TT 
Î u(x, + rcos8, y + rsin6)d@ -| u(x + ECOS 6 , y, + esin 0)d06. 


la frontière de U. 


3. Les fonctions harmoniques dans U possèdent d'autres pro- 
priétés. En particulier, toute fonction harmonique dans U est 
nécessairement indéfiniment dérivable (et ses dérivées sont 
bien sûr harmoniques) et analytique dans U. 


4.3 L'équation de Laplace et le problème 
de Dirichlet dans un cercle 


On se propose de résoudre le problème suivant : trouver les fonc- 
tions harmoniques U (c'est-à-dire solutions de l'équation de 
Laplace) dans le disque de centre O de rayon R qui se prolongent 
par continuité sur la circonférence de centre O et de rayon R et 
vérifient : 


U(Rcos0 ,Rsin0 ) = f(8) 


où fest une fonction continue 2r-périodique. 


Autrement dit, on cherche les fonctions harmoniques U prenant 
des valeurs données sur la circonférence de centre O et de rayon A. 


4.3.1 Unicité d’une éventuelle solution 
par le principe du maximum 


Supposons que l'on dispose de deux solutions U; et U, de ce 
problème de Dirichlet. Celles-ci sont harmoniques dans le disque de 
centre O et de rayon R, se prolongent par continuité à la circonfé- 
rence de centre O et de rayon R où elles sont égales. 


D'après la proposition 3, elles sont égales. 
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4.3.2 Existence d’une solution par la méthode 
de séparation des variables 


Puisque l’on travaille sur un cercle, nous allons utiliser des coor- 
données polaires de centre O. Au lieu de chercher directement une 
fonction harmonique (x, y) — U(x, y) solution du problème, nous 
allons poser 


U(x, y) = U(r cosä rsin0 ) = u(r, 0) 


et déterminer d’abord u. 


Vu l'expression du laplacien en coordonnées polaires, on sait que 
U vérifie AU = 0 si, et seulement si, u vérifie la relation : 


d2u 1 
—— + — 
: 


ar? 


ou, 1 du | 
dr 1 d®@ 

De plus, cette fonction u devra vérifier la condition au bord 
u(R,8)=#f(0). 


# On commence, comme d'habitude, par chercher les solutions 
sous forme multiplicative 


(5 80)—u(rv(8 ), 
ce qui conduit à résoudre deux équations de la forme : 
r2u’{(r)+ru’(n+2Au(rn = 0 ; 
v’’(8)}-Av(8)=0 


où À désigne un nombre réel. 


Pour tenir compte du fait que la fonction obtenue doit être de 
classe C? sur le disque de centre O et de rayon P, il faut bien sûr 
que : 


v(0) = v(2n), 
v’(0) = v’(2x), 
v’’(0) = v’’(27) 


(ce qui conduit encore à un problème de Sturm-Liouville). Ces 
dernières conditions ne peuvent être réalisées que si, et seulement 
si, À est de la forme - n? où n désigne un nombre entier naturel non 
nul. 


L’équation en u est une équation d’Euler, dont on cherche des 
solutions dépendantes sur 10, + + [ sous la forme r- r4, de sorte 
que les solutions obtenues pour les deux équations précédentes 
sont, avec r> 0: 


v (8) = A, cos(n0 ) + B, sin(n0) ; 
utrieC. ri Dr 


Pour que u se prolonge par continuité sur [0, R], il faut de plus 
que D, = 0. 


# L'équation étant linéaire, on cherche donc, de façon formelle, des 
solutions de la forme suivante : 


u(r, 0 ) = À; DE Ja, cos(n0 )+ B, sin(n0 )). 


Pour que la condition au bord du cercle soit vérifiée, il faut enfin 
que, pour r=R: 


f(0) = Dr cos(n6 )+ B, sin(n6 )). 


On choisira donc pour nombres A, et B, les coefficients de 
Fourier de fet l’on est finalement conduit à envisager pour solution 
de notre problème la fonction : 


+ co 


1 27 r\9/1 27 
u(r, 0 ) = [L f (6) dt+ >. (+) GL, f (t)cosnt dt cosn0@ 
n=1 
1 2x à 
| f () sinnt dt sin(n0 )). 
0 


Cette expression peut également s'écrire : 


u(r, 0 ) = L [rc + F} ()'eosnce 0 ar 


n=1 
Tr D Car 
= 1F f(0-t) 5 cost ét : 


En effet, la fonction f étant continue et 2x-périodique est bornée, 
et il est clair que, pour 0<r<R, la série figurant ci-dessus sous le 
signe intégral est uniformément convergente sur [0, 2x], ce qui 
justifie la permutation des symboles de sommation et d'intégration. 
La seconde expression obtenue résulte alors de la première à l'aide 
du changement de variable t” = 0-t. 


# On s'intéresse maintenant à la fonction P suivante, définie pour 
tout couple (x, t) de nombres réels tels que — 1 < x< 1 et appelée 
noyau de Poisson du problème de Dirichlet relatif au cercle : 


+ co 


1,3 D x? cosnt. 


P(X, t) = Bet 


n=1 


On remarque d’abord que cette fonction est paire en la variable t 
et que son intégrale sur [0, 2x] est égale à 1 (on peut en effet intégrer 
terme à terme en t par convergence uniforme de la série). De plus, 
on a l'expression plus simple : 


+ 00 


+ co 
1 4 n np Hyn 
P(X, t) = set 2% cosnt = 3:+7 Re > Greit) 
n=1 hz=1 . 
it 2 
= TR LA à 1 X 


2x %  1-xet 2n(x2-2x cost+]1) 


Finalement, on dispose des deux formes suivantes pour u (r, 6): 


27 
ur. 0 ) = [ f(0-t) P(R rat 


Pet tm dt 
7 Jo os ETAT: cost+R?2) 


u(r, 0 ) = ['r P (& e-t}at 


“ft _ d 
= dt, 
[ D > -2rR cos(0—t)+R 2) 
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Revenant en coordonnées cartésiennes, cette formule s'écrit 
aussi : 


R2-x2- y2 
2r(xX2+y2-2R (x cost+ y sint)+R2) 


27 
U(Xx, y) T [L f(t) 


Il reste maintenant à voir si cette fonction est bien solution du 
problème de Dirichlet que l’on s'était donné au départ. 


e Cette fonction est continue, de classe C?et peut se dériver sous 
le signe intégral pour : 


0O<r= /x2+y2<R 


d'après les théorèmes de continuité et de dérivation sous le signe 
intégral. En effectuant le calcul, on vérifie, de plus, facilement que : 


R2-x2- y2 


27 
AU(X, y) = ft ER EL RTE TO TETE 
alé [ ‘0 2r(x2+ y2-2R (x cost+ y sint)+ AR?) 


Jat = 0 


e Cette fonction se prolonge par continuité sur la circonférence 
de centre O et de rayon P, et u{(r, 8) tend bien vers f(6,) lorsque 
(r, 6) tend vers (PR, 65). En effet, puisque l'intégrale du noyau de 
Poisson vaut 1,ona: 


R2-r2 


— dt. 
2r(r2-2rR cost+R 2) 


27 
u(r, 0 )—-f (60) = f (A(8- t)- (@5)) 


La fonction f étant continue et 2rx-périodique est évidemment 
bornée sur R (par un nombre M) et uniformément continue sur R. 


Pour tout nombre € > 0, il existe donc un nombre «> 0 tel que : 


(x-y<20)=(FCO-F(YIE e). 


Supposons donc [8-6 <aetR sina<r<R et montrons que 
u(r, 0) -f(0 ) est alors inférieur en valeur absolue à 3e. 


A cet effet, on découpe l'intervalle d'intégration [0, 2x] en trois 
intervalles, allant de 0 à a, de œ à 2n- «et de 2x -a à 2n. On a 
clairement : 


If(8-D-F(8) <E 


sur le premier et le dernier de ces trois intervalles et, le noyau de 
Poisson étant positif et d'intégrale égale à 1 sur [0, 2x], on peut donc 
majorer par £ les intégrales sur ces premier et dernier intervalles, de 
sorte que l'on a: 


[u (r0)-f (65) 


R2-r2 


27x- a 
<2 ROSÉ ON M 
e+[ ESP Te 26m costs À 9 


< 2€+ 


2M(R2-r2)[{27-a dt 
2T « (r-Rcost)}2+R 2sin ?t 


Maintenant, il est clair que le dénominateur dans cette dernière 
intégrale est supérieur ou égal à R?sin?œ (c'est évident pour 
a<t<n/2 et 3n/2<t<2r-a en minorant par 0 le premier 
terme situé au dénominateur ; c'est également vrai pour 
x/2<t<3x/2 car r- Rcost>=r>z Rsinæ.Finalement, on voit que, 
quitte à choisir rencore plus proche de P, on a bien le résultat voulu, 
c'est-à-dire : 


[u(r, 8 )-f (8) <2e+ 


Le problème de Dirichlet est ainsi résolu sur le cercle. 
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On laisse enfin au lecteur le soin de vérifier qu'ici encore la solu- 
tion u obtenue pour le problème de Dirichlet sur le cercle dépend 
continüment de la donnée initiale fet de tirer les conséquences de 
ce résultat. 


Principe de la moyenne et harmonicité sur un ouvert. | 


On est maintenant en mesure de prouver (en restant dans le 
cas de deux variables, mais la démonstration se généralise sans 
grande modification au cas général) que toute fonction u conti- 
nue sur un ouvert U et y vérifiant le principe de la moyenne est 
harmonique (réciproque de la proposition 4). 


Considérons un point M, appartenant à U. Comme U est 
ouvert, il existe un nombre réel R > 0 tel que la boule fermée de 
centre M, et de rayon R est incluse dans U. Introduisons alors, 
sur cette boule fermée, la fonction v définie par : 


v(X Y) = 


R2 - x2 - y? 


T 
u(X) + Rcost, RE 2 Route sono 0e 
[L (Xo + Yo + rca y?-2R(xcost+ ysint)+R?) | 


(On a repris la formule précédente en changeant f{t) en 
u (x + R cost, yo + R sin). 


D'après ce qui précède, cette fonction v est continue sur la 
boule fermée de centre M, et de rayon R ; elle coïncide avec u 
sur la circonférence de centre M, et de rayon R et elle est har- 
monique dans le disque de centre M, et de rayon PR, donc elle y 
vérifie le principe de la moyenne d'après la proposition 4. 


Les fonctions u et v sont donc égales sur la circonférence de 
centre M, et de rayon R et vérifient le principe de la moyenne 
dans le disque de centre M, et de rayon R, donc le principe du 
maximum (cf. encadré du paragraphe 4.2). Il résulte alors de la 
proposition 3 que u = v sur la boule fermée de centre M, et de 
rayon R et u, comme v, est donc harmonique dans le disque de 
centre M, et de rayon R, et par conséquent dans U puisque Mo 
est un point quelconque de U. 


4.4 L'équation de Laplace et le problème 
de Neumann dans un cercle 


4.4.1 Généralités 


On se propose de résoudre le problème suivant : trouver les fonc- 
tions harmoniques U (c'est-à-dire solutions de l'équation de 
Laplace) dans le disque de centre O de rayon R qui se prolongent 
par continuité sur la circonférence de centre O et de rayon À et y 
admettent en tout point M (AR cos0, R sin@) un vecteur dérivé suivant 


A ., 
le vecteur unitaire normal n (cos@ ,sin@ ) vérifiant la condition 
dite de Neumann : 


OU (RcosO ,Rsin0) = (0) 
on 


où fest une fonction continue 2r-périodique. 


Autrement dit, on cherche les fonctions harmoniques U dont la 
dérivée normale suivant la circonférence de centre O et de rayon R 
est donnée. 
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Remarque préliminaire. Rappelons que le flux du gradient 
d'une fonction harmonique au travers d’un contour (ou d'une sur- 
face) fermé est nul. En effet, la première formule de Green (8 4.2) 
s'écrit : 


+ —7 av 
J] Vu -Vv dxays [] u Av dxdy = u— ds. 
D D on 


oD+ 


Appliquée avec u = 1 et v= U, où U désigne une solution du 
problème de Neumann envisagé ci-dessus, on en déduit que la 
fonction fvérifie nécessairement la condition : 


2 27 
$ CRU = Ï *JU(Rcose , RsinO )R d0= a f(0 )d@ = 0. 
0 ün 0 


On verra que cette condition nécessaire d'existence des solu- 
tions est, en fait, suffisante. 


4.4.2 Unicité à une constante additive 
près d’une éventuelle solution 


Supposons que l’on dispose de deux solutions U; et U, de ce 
problème de Neumann. Celles-ci sont harmoniques dans le disque 
de centre O et de rayon R, se prolongent par continuité à la circonfé- 
rence de centre O de rayon R où elles vérifient pour tout 0: 


CAUA : aU, : 
— (Rcos0, Rsin0 ) = — (RcosO , Rsin0 ). 
on on 


D'après la formule de Green que l’on vient de rappeler, appliquée 
avec u = V= U; - U,, on a: 


—2 
(ll V (U,-U,)(x ydxdy = 0. 
B(0,R) 


Comme U; - U, est de classe C2, son gradient est nul sur la boule 
de centre O et de rayon R, et U, - U, y est constante. 


Réciproquement, on vérifie que si U est solution du problème de 
Neumann précédent, alors U + Cte l'est également. 


4.4.3 Existence d’une solution par la méthode 
de séparation des variables 


Puisque l’on travaille sur un cercle, on utilise des coordonnées 
polaires de centre O. Au lieu de chercher directement une fonction 
harmonique (x, y) — U(x, y) solution du problème, on pose : 


U(x, y) = U(r cos, rsin9 ) = u(r, 0 ). 


et on détermine d’abord u. 


& Les calculs conduits au paragraphe 4.3.2 amènent encore à 
considérer (au moins formellement) à une constante additive près : 


u(r, O ) = 25 cos((n0 + B, sin(n0 ))). 


Tenant compte de la relation immédiate : 


aU 
— 


(r cos8 , rsin0 ) = ou (r,0) 
“a or 


on voit qu'il faut enfin, pour que la condition de Neumann soit véri- 
fiée, que : 


(8) => F (A, cos(n6 )+ B, sin(n6 )). 


On choisira donc pour nombres nA,/R et nB,/R les coefficients 
de Fourier de f(on a A) = 0 à cause de la condition nécessaire sur f 
obtenue initialement) et l'on est finalement conduit à envisager 
pour solution de notre problème la fonction : 


u(r,0)= 

S° i 27 27 & . 

>, (r) (A f (t)cosntdt cosn0+ R f(bsinnt dt sin n6 ). 
n= 1 R nr 0 nT b 


Cette expression précédente peut également s'écrire : 


au 8)= [rc » HE cosn(6- t}dt 
" n=1 


T r® n 
& L f(8-t) ÿ (+) À cosnt dt. 
n=1 


En effet, la fonction fétant continue et 2r-périodique est bornée, 
et il est clair que, pour 0<r<R, la série figurant ci-dessus sous le 
signe intégral est uniformément convergente sur [0, 2x], ce qui 
justifie la permutation des symboles de sommation et d'intégration. 
La seconde expression obtenue résulte alors de la première à l'aide 
du changement de variable t”=0-t. 


On s'intéresse maintenant à la fonction P suivante, définie pour 
tout couple (x, t) de nombres réels tels que — 1 < x< 1 et appelée 
noyau de Poisson du problème de Neumann relatif au cercle : 


+ co 
n 
FxDS 7% _ cosnt . 


n=1 


Par dérivation terme à terme de cette série entière en x, on a 
donc : 


+ 00 as 
5x6 t) = L y. x-1 cosnt = ! Re Le" > eee y 


n=1 n=0 


eit cost- x 


de 
rm 1-xelt n(x2-2x cost+1) 
Par intégration, on obtient alors : 

LL es 

P(CXx, t) = 5x In(x£ —-2x cost+1). 
Finalement, on dispose de la formule suivante pour u (r, 6): 
R|"#(0-0 PL, td 
0 cd. FR ) : 


T VER 2 
2 | f (6-0 in(£ 2rR cost+R at 
27 0 


u(r,0 ) 


R?2 


ou, en coupant le logarithme et en utilisant la nullité, l'intégrale de f 
sur [0, 2x] : 


T 
u(r, 0) = - 2 | f(8-t) In(r2-2rR cost+R?)dt 


u(r, 0) = = À [re In(r2-2rR cos(0-t)+R?2)dt. 
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Revenant en coordonnées cartésiennes, cette dernière formule 
s'écrit aussi : 


T : 
U(x, y) = - s [ f(t)In(x2+y2-2R (xcost+ ysint)+R2)dt. 


Comme pour le problème de Dirichlet, il resterait maintenant à 
voir si cette fonction est bien solution du problème de Neumann 
que l’on s'était donné au départ. 


4.5 Le potentiel newtonien 
et l'équation de Poisson 


On sait, en mécanique, que les forces de gravitation exercées sur 
un point quelconque par un corps de volume V'et de masse volu- 
mique p proviennent d'un champ de gravitation g dérivant d'un 
potentiel de gravitation &. En un point Molxo, Yo, 2) de l'espace, le 
potentiel de gravitation exercé par ce corps D est égal à: 


eu efff, 


- G 


p (x y, z)dxdydz 


(X- Xp)? +y- Yo)? + (a 25) 


p dV 


où G désigne la constante de la gravitation (égale à 6,67.107!1 S.1.) et 
où la fonction de masse volumique P évidemment nulle à l’exté- 
rieur du Corps considéré, sera supposée continue sur celui-ci et de 
classe C! à l’intérieur de celui-ci. 


On établit alors en mécanique les relations suivantes : 


— 
g=-Vo, 
div g=-4rGp. 


Puisque la divergence du gradient d’une fonction numérique de 
plusieurs variables n'est autre que son laplacien, on déduit de ces 
deux relations que : 


AD = 4nGp. 


On dit alors que © est solution de l'équation de Poisson et 
puisque p = 0 à l'extérieur du corps considéré, on a A®= 0 en tout 
point extérieur de celui-ci, ce qui signifie que le potentiel de gravita- 
tion ® créé par un corps est harmonique à l'extérieur de celui-ci. 


Exemple : potentiel newtonien créé par une boule lors- 
que p ne dépend que de r. 


On considère la boule B de centre O et de rayon R et l'on suppose 
que sa masse volumique p{x, y, z) = pr) ne dépend que de 


r= Nx2+y2+7z2. 


Pour évaluer le potentiel newtonien qu'elle crée en un point M, 
nous allons supposer (ce qui ne change rien vu la symétrie du pro- 
 blème) que celui-ci appartient à l'axe Oz et qu'il a donc pour coordon- 
nées (0, O0, r). On calcule alors l'intégrale précédente, mais en passant 
en coordonnées sphériques, ce qui revient à poser : 


x = rsin0 cos , 
y = rsin@ sin® , 
Z = f C08Û . 


Le jacobien de cette transformation étant égal & à r? sin6, on obtient 
donc : 


p (ndxdydz 


® (Mo) = -G ———}"—@ _—_—_—_———— 
Là IT. Jx2+y2+(2-10)? 


aff * r?p (Dsin0 dé de dr 


x ,fr2-2rr cos0+ rà 


On intègre d'abord en #, puis un simple calcul conduit à : 


® (M5) = -2s[0 POI Jr2-2rrcos8+r2 ]6 dr 


= ef rp(Ilr+ ro -|r-"lidr. 


Il convient alors de distinguer suivant que M, est extérieur ou inté- 
rieur à la boule B. 

e Si M, est extérieur à la boule, c'est-à-dire si Rr ZAR, on a en 
notant m la masse de la boule : 


ArG 


® (Mo) = - à e p(n) dr= - 09, 


On remarque que ce potentiel est tu à celui que créerait la même 
masse m située ponctuellement au centre 0 de la boule B considérée. 


e Si Mb et intérieur à la boule, c'est-à-dire si nr <R,0on a: 
® (Mo) = - + la "Ê p(r) dr- axe [° rp(r) dr. 
0 


On vérifie que © est continue sur l'espace, ÿ compris sur la sphère 
de centre O et de rayon R et l'on vérifie par un simple calcul de déri- 
vées que : 

— AD = 0 à l'extérieur à la boule (@ est harmonique à l'extérieur de 
la boule), 

— AD = 4nGp à l'intérieur de la boule, et en fait partout puisque 

. p = 0 à l'extérieur de la boule (® vérifie donc l'équation de Poisson). 


orie spectrale 
saration des 


Pour de nombreuses équations aux dérivées partielles abordées 
précédemment, c'est la méthode de séparation des variables 
(consistant à rechercher d'éventuelles solutions de la forme : 


(X y) uCv(y), 


puis à en former des combinaisons linéaires éventuellement infinies 
en un sens à préciser) qui a souvent permis de déterminer les solu- 
tions recherchées. 


Le succès théorique (et également pratique pour de nombreux 
cas particuliers simples) de cette méthode repose sur les fonde- 
ments que nous exposons maintenant. On a vu, dans les paragra- 
phes précédents, que la recherche de solutions de la forme 


(x y) = u(v(y) 


vérifiant des conditions aux limites imposées conduisait en général 
à la résolution d'équations différentielles, l’une d'inconnue u, l'autre 
d'inconnue v, l’une de ces deux fonctions au moins étant de plus 
astreinte à vérifier des conditions aux limites. Comme on l'a déjà 
indiqué, il s’agit en fait de résoudre ce qu'il est convenu d'appeler 
un problème de Sturm-Liouville, que nous définissons ci-dessous 
de façon plus précise. 
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Définition 2. 

Soient p, q, rtrois fonctions à valeurs réelles continues sur un 
segment [a, bl], la fonction p étant de classe C' sur La, bl et les 
fonctions p et r étant à valeurs strictement positives sur [a, b]. 

Soient (A, B) et (C, D) deux couples de nombres réels dis- 
tincts de (0,0). 

On appelle alors problème de Sturm-Liouville régulier associé 
à ces données la recherche des nombres réels 1 et des solutions 
non nulles y de l'équation différentielle du second ordre : 


(p{t)y') +q{t)y+2Ar(t)y=0 
qui vérifient de plus les conditions 
Ay (a) + By'(a) = 


et 
Cy (b) + Dy'{(b) = 


Exemple : Le plus simple des problèmes de Sturm-Liouville déjà 
abordés est le suivant : 
y’’+y = 0 


. avec y (0) = 0 et y (x) = 
dont les solutions sont les nombres 


= n? 
et les fonctions associées 
y(t) = C, sinnt 


où n désigne un nombre entier naturel non nul et C, une constante 
réelle. 

De plus, toute fonction f, continue et de classe Gé" par morceaux sur 

. [0, x] qui vérifie de plus les deux conditions aux limites f(0) = f(x) = 
peut être développée en série normalement convergente sous la 
-: forme : 
+ © 
2 
f(x) = . br 


n=1 


ET sin(nt) dt sinnx. 
0 


Cela n'est ici qu'une conséquence du théorème de convergence 
normale de Dirichlet sur les séries de Fourier, appliqué à la fonction f 
supposée prolongée sur R de façon impaire et 2x-périodique. 


Nous verrons plus loin que les résultats de ce problème parti- 
culier se généralisent en fait aux problèmes de Sturm-Liouville 
réguliers. 


Auparavant, donnons une interprétation de notre problème de 
Sturm-Liouville en termes de valeurs propres et de fonctions 
propres. À cet effet, on désignera par : 


+ S l’espace vectoriel des fonctions de classe C? de [a, b] dans 
R vérifiant de plus les deux conditions aux limites : 


Ay (a) + By'(a) = 0 


et 
Cy (b) + Dy"{b) = 


+ L l'application linéaire définie de l’espace vectoriel S défini ci- 
dessus dans l'espace vectoriel C([a, b], R) des fonctions conti- 
nues de [a, b] dans R par la formule : 


Ly (D = [(PDY (0) + gay]. 


EG 


Le problème est alors de rechercher les nombres réels À pour 
lesquels existent des solutions non nulles y appartenant à S de 
l'équation Ly = Ày. 
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Si l’on met de côté le fait que L n'est pas un endomorphisme 
d'espaces vectoriels, mais une application linéaire de S dans 
C([a, b], R}), on voit qu'il s’agit de rechercher les valeurs propres À 
et les fonctions propres associées y de L. 


Proposition 5. 


Dans le contexte introduit dans ce paragraphe 5, les sous-espaces 
propres de l'application linéaire L sont tous de dimension 1. 


Démonstration. © 


Un sous-espace propre E (À) est par définition de dimension au 
moins 1 et il est ici de dimension au plus 2, car l'égalité : 


Ly= Ày 
implique que y est solution de l'équation différentielle : 


py"+p'y'+{q+2r)y=0, 


et l’on sait par le théorème de Cauchy-Lipschitz que les solutions de 
cette équation différentielle homogène du second ordre forment un 
espace vectoriel de dimension 2. 


Enfin, la dimension d’un sous-espace propre E (À) est exactement 
égale à 1 ici car, si elle était égale à 2, les éléments de E (1) seraient 
exactement (par égalité des dimensions) les solutions de l'équation 
différentielle 

py"+p'y'+{(q+r)y=0; 
autrement dit, toute solution de cette équation appartiendrait à S et 
vérifierait les deux conditions aux limites : 
Ay (a) + By”(a) = 
Cy (b) + Dy'(b) = 


Cela est impossible car la solution vérifiant les conditions initiales 


y (a) = 
et 
y (a) = 
ne vérifie pas la première relation puisque (A, B)Z0. () 


Proposition 6. 


Dans le contexte introduit dans ce paragraphe 5, l'application 
linéaire L est autoadjointe pour le produit scalaire défini sur S par: 


<fg>= [rc g(®) rt) dt. 


Démonstration. © 
L'énoncé signifie simplement que, pour tout couple (7, g) de S, on 
a l'égalité : 
slide <f,l9s, 


Celle-ci se vérifie par intégration par parties. On a en effet : 
b 
<Lfg>= -| [(pCOF "(0)" + (0) Og(bdtl 
b b 
= [ p(0 F(0g a+ | POr' bg (bdt- | q(Of (Og(Odt 
a a 


et, si l’on suppose les constantes À et C non nulles (on étudie plus 
simplement encore le cas où À ou Cest nulle), il vient : 


<Lfg>= D C p(b) fa(b)- p(a)fg(a) 


b 
+ [ec ‘g'(t - q(tfg(tidt.. 
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On voit maintenant que f et g jouent des rôles symétriques et, 
en permutant fet g puis en remontant les calculs, on a bien l'égalité 
voulue. 


Ce résultat est intéressant car on connaît déjà, au moins dans le 
contexte des espaces euclidiens (il s’agit des espaces de dimension 
finie sur R munis d’un produit scalaire), le théorème spectral clas- 
sique concernant les opérateurs autoadjoints : 


— les valeurs propres sont réelles ; 

— les sous-espaces propres sont deux à deux orthogonaux ; 

— la somme directe des sous-espaces propres est égale à 
l'espace entier. 


On démontre ici de la même façon les deux premiers points, de 
façon élémentaire, mais le troisième point doit, par contre, être 
précisé, et sa démonstration est beaucoup plus difficile : c'est le 
théorème spectral d'Hilbert-Schmidt, dont on donne l'énoncé. 


Proposition 7. 


Dans le contexte introduit dans ce paragraphe 5, les valeurs 
propres de l'application linéaire L définies sont simples et peuvent 
être rangées en une suite croissante 


A0 <<<: <À,<: 


tendant vers + ce. 


Si l'on désigne par e, une fonction propre unitaire (autrement dit 
de norme 1 au sens de la norme dérivant du produit scalaire précé- 
dent) appartentant à S et associée à la valeur propre À, de L, la suite 
(e,) est orthonormale et l’on a alors pour toute fonction f apparte- 
nant à S (et, plus généralement, pour toute fonction continue et de 
classe C' par morceaux sur [a, b] vérifiant les deux conditions aux 
limites : 


Af(a) + Bf'(a)=0 


et 
Cf{(b) + Df'(b) = 0 


l'égalité suivante, où a< x< b : 


+ co 


HD = Y<e,sfs 6,0%. 


n=0 


De plus, la convergence de cette série est absolue et uniforme sur 
[a, b]. 


Dans les exemples étudiés dans les paragraphes précédents 
(8 2 à 8 4), c'est en fait ce théorème que l’on a démontré élémentai- 
rement à chaque fois, étant entendu que celui-ci est susceptible de 
généralisations (théorème dit de Weyl-Stone-Titchmarch-Kodaira) 
permettant l'étude des problèmes singuliers de Sturm-Liouville qui 
apparaissent lorsque le segment d'étude [a, b] est remplacé par un 
intervalle réel quelconque 1, non nécessairement fermé ou non 
nécessairement borné. 


Les suites des fonctions (e, ) qui s’introduisent ainsi conduisent à 
ce qu'il est convenu d'appeler en mathématiques les fonctions 
spéciales, et qui comprennent notamment toutes les familles de 
polynômes orthogonaux usuelles (Legendre, Hermite, Laguerre, 
Tchebichev, etc.) et les fonctions de Bessel qui s'introduisent avec le 
problème de Sturm-Liouville singulier associé à l'équation suivante 
(où v désigne un nombre réel arbitrairement donné) : 


v 2 
(xy’ y - y+Àxy = 0. 


On donne maintenant, pour finir, l'exemple d'un problème de 
Sturm-Liouville singulier conduisant à l'utilisation de ces fonctions 
de Bessel. 
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Exemple : membrane vibrante d’un tambour circulaire 


M Considérons la vibration de la membrane d'un tambour circulaire 
que nous assimilerons au repos au disque de centre O et de rayon À du 
plan. La cote U{x, y, t) de cette membrane au point de coordonnées 
(x, y) à l'instant t vérifie l'équation : 


a2U di _ 
= V + 


2 ox ay 


— avec la condition aux limites : 
U(Rcos® , Rsin6,t)=0 


pour 8 quelconque et = 0, 
— avec les conditions initiales : 


U(Xx y,0) = f(x) 
et 


oU 
Er (X y, 0) = 0 


où la fonction fest supposée suffisamment régulière et nulle sur la cir- 


conférence de centre O de rayon A. 
Nous allons rechercher les solutions en passant en coordonnées 


: polaires, en prenant la nouvelle fonction inconnue définie par : 


u(r, 8, t) = U(rcos8, rsin8 ,t). 


Celle-ci vérifie donc, outre la condition aux limites u (A, 8, t) = 0 et 


- des conditions initiales faciles à préciser, l'équation aux dérivées par- 
 tielles suivante : 


du _ 2 = 


1 d2u 
EYA FT ) 


L'on, 1 eu 
r dr r2 202 


Pour simplifier, nous en chercherons seulement les solutions possé- 


dant la symétrie de révolution : le dernier terme de l'équation ci-dessus 


est alors identiquement nul et on applique dans ce contexte la 


méthode de séparation des variables. 
On recherche donc des solutions de la forme : 


(5 D — un) v(P) 


et telles que : 
u(R)= 0. 


ILest facile de voir qu'une telle fonction est solution si, et seulement 


si, uet v vérifient les équations : 


ru ”’+u’+ Àru = 0 
V” = JV 


avec les conditions u =- Av? et u(R)= 0. 


On voit facilement que la fonction u s'obtient par résolution d'un pro- 
blème de Sturm-Liouville singulier sur le segment [0, AR] puisque la 
fonction rs'annule en 0. 

La condition aux limites est simplement u (R) = 0 (la condition en O 
se limitant à exiger la continuité de u en ce point singulier pour le pro- 
blème). 


M Ce problème de Sturm-Liouville est en fait bien connu et on peut le 
résoudre explicitement. On vérifie tout d'abord que le nombre réel À 
est nécessairement positif. En multipliant la première de ces équations 
différentielles par u, on obtient en effet en intégrant sur [0, R]: 


[eu (nu”’(n+ u(n) u’ (n+ ru 2(n))dr = 0. 
0 
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Par intégration par parties sur les deux premiers termes, il vient 
alors : 


D LE mi 
[ru (Du (nlo - … u'£ (ndr+2 ni (dr = 0. 


Et comme le crochet est nul, la positivité de À en résulte. 


On posera donc À = w? et on constate que la fonction r— u(r/@) 
est solution de l'équation de Bessel : 


y"+y"+ry=0. 


Les seules solutions continues en 0 de cette équation de Bessel 
d'indice 0 sont de la forme 


r— CJ(n) 


où C'est une constante réelle arbitraire. Ainsi, on a : 
u(r) = CJ(@r) 


et comme u doit s'annuler en r = R, on voit que les nombres œR sont 
des zéros de J. Comme on sait que les zéros des fonctions de Bessel 
(et de J en particulier) forment une suite, que nous notons ici (x,), de 
nombres réels strictement positifs, on voit que æ est nécessairement 
de la forme x,/R. 


D Afave EN 
Références 


Par conséquent, les solutions multiplicatives de l'équation 
(r D u(rn v(t) 
vérifiant u (A) = 0 sont de la forme 
(0 = Jiéfrx,/RYCA, cos(vtx, /R) +B, sin(vtx, /R)). 


# Pour achever le problème, il reste à chercher des solutions par 
superposition, donc de la forme : 


rx vtx} _ vEx; 
u(r, t) = Z0( (ar cos” R +B, sin A) 


et vérifiant les conditions initiales requises. 


On choisit à cet effet B, = 0 et il reste à savoir si l’on peut déter- 
miner des nombres A, tels que l'on ait pour 0<r<R : 


60 = SA d( ): 


Un tel choix des A, est effectivement possible lorsque Ÿ est de 
classe C? sur [0, R ] (maisil suffirait en fait qu’elle soit continue et de 
classe C | par morceaux sur [0, R ]) et nulle en R, d'après l'extension 
signalée plus haut du théorème de Hilbert-Schmidt au cas de notre 
problème singulier de Sturm-Liouville. 
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es systèmes dynamiques n'ont été étudiés en tant que tels que assez tardi- 
. vement. Ils sont néanmoins apparus assez tôt dans l’histoire scientifique, 
puisque l'on peut les reconnaître dans les premiers travaux de la mécanique 
donnant lieu à des équations différentielles. 

Schématiquement, un tel système est la donnée d’une loi d'évolution qui, à 
partir de conditions initiales, détermine le futur d’un phénomène. Le paradigme 
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en est l'équation différentielle, qui exprime une loi régissant, elle-même, l'évolu- 
tion temporelle d'un phénomène convenablement paramétré. Cette loi déter- 
mine l’évolution du système lorsque les paramètres sont connus à un certain 
instant. Sous cette forme, le système dynamique ne peut rendre compte que 
d’une loi déterministe. | 


La résolution explicite, ou même approchée, d'une équation différentielle, est 
en général impossible. Dans une large mesure, l'étude des systèmes qui nous 
occupent vise à formuler les termes d’une étude qualitative des phénomènes. 


Dans le cadre de cet article, nous nous contenterons d'introduire le langage 
_ nécessaire, en restant dans un cadre élémentaire. Nous n'aborderons pas vrai- 
ment le problème fondamental de la perturbation dans un système dynamique 
ou d’un système dynamique, qui fait l’objet d’un article séparé, dont la lecture 
suppose la connaissance des notions développées dans le présent article 


Le but de la partie 1 est de fournir des exemples d'équations différentielles sur 
lesquels est présenté informellement le langage. Celui-ci sera précisé dans le 
cadre de la partie 2, où sont aussi fournis les outils élémentaires sur les équa- 
tions différentielles. La partie 3 donne l’occasion de mettre en œuvre ces outils. 
La partie 4 dégage, sans y insister, les propriétés fondamentales, qui permettent 
d'introduire les systèmes dynamiques discrets. Dans la partie 5, on trouvera trai- 


tés des exemples très simples de systèmes dynamiques dans R. 


1.1 Exemple d’une équation vectorielle 
linéaire à coefficients constants 


Soit À une matrice carrée réelle, de taille p. On considère l'équa- 
tion différentielle : 


(E) X7 = A 


Les solutions de cette équation différentielle sont, par définition, 
les fonctions X de R vers me, 1(R), satisfaisant à l'égalité : 


VteR X'(t) = AX(t). 


Nous savons parfaitement résoudre cette équation différentielle 
linéaire à coefficients constants, de taille p, d'ordre 1. La solution 
générale en est de la forme : 


te X(H = etAC, 


où Cest un élément arbitraire de M, ,(R). 


æ Si l’on cherche plutôt à résoudre le problème de Cauchy où 
l’on impose en outre à X de prendre, pour la valeur t de la variable, 
la valeur X,, on voit que C est déterminé par l'égalité : 


X, = eAC, 
soit C = ekAX, 
et donc: 
VteR X(t) = ett-Ax,, 


On remarque au passage qu'il existe une, et une seule, solution au 
problème de Cauchy. 


Nous prenons désormais f, = 0. La variable t désignant souvent 
une variable temporelle, nous dirons que la condition initiale est 
fixée à l'instant 0. Ainsi : 


VtER X(t) = efAX,. 


L'application tr efAX, est à valeurs dans me, 1(R) . Pour se con- 
former à l'usage, on identifiera librement cet espace vectoriel à RP. 
Ainsi, t+> efAX, est un arc paramétré de RP, appelé trajectoire de 
Xo- Il s’agit donc d'une notion cinématique, puisque cette applica- 
tion décrit la position de X à l'instant f. 


Il arrive que l’on s'intéresse seulement à l'ensemble des positions 
atteintes par X. On appellera orbite de X, cette courbe, qui n'est rien 
d’autre que l’image de l'application X. 


En réalité, X (t) dépend, outre de t, de la position initiale X$4. Pour 
bien marquer cette dépendance, nous poserons : 


Ainsi, t- ÿ(t, X,) est la solution du problème de Cauchy : 
X°= AX : X(0) = X;. 


Pour rendre plus parlante l'application #, on peut considérer que 
l'on soumet RP, constitué de points matériels, à une loi reflétée par 
l'équation différentielle X”7 = AX. Si l’on observe un point X, de 
RP, c'est-à-dire un point matériel qui, à l'instant initial, occupe cette 
position, sa position à l'instant f sera @(t, X,). En particulier, si l'on 
se donne un certain nombre de valeurs de XQ («position du 
système » à l'instant initial), la famille des ÿ(t, Xç) pour un certain t 
décrit la position du système à l'instant t, c'est-à-dire, en quelque 
sorte, sa photographie à cet instant. Si l’on fixe une seule valeur de 
X, et que l'on fait au contraire varier t, on obtient l'évolution tempo- 
relle d’un point matériel. 


L'application 4, qui contient simultanément l'ensemble des infor- 
mations précédentes, est appelée flot de l'équation différentielle 
X' = AX. 
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x Pour proposer quelques illustrations, nous fixons désormais À, 


en prenant p=2et À = | $ : | On obtient, de suite : 


bCE Xo) = | . 0 pe 


e2t 


Posons : 
X = (X1: X2) et LXGA Xo) = (d1CE Xo) CE Xo)): 


Nous avons ici identifié les matrices colonnes de taille 2 aux élé- 
ments de R2. Ainsi : 


d: = ex, ; hp = e?tx,. 


L'étude de cet arc paramétré ne présente aucune difficulté : 
ix>0 et x, >0, l'orbite est l'arc de la parabole ÿ, = =262 
— si x,>0 et x,>0, l'orbite est l'arc de la parabole ÿ, = 29 
1 


limité à la portion de cette parabole incluse dans le premier qua- 
drant (ouvert) ; 

— on obtient trois autres études identiques selon les signes de x; 
et x, supposés non nuls ; 

— si x ou x est nul, sans que les deux le soient, on obtient une 
demi-droite ; 

— enfin, si, x = xX = 0, il s’agit du point (0,0). 


æ Visualiser globalement l'application Ÿ présente quelques diffi- 
cultés. 


e Une première idée est de visualiser l'orbite d'un point X ; mais 
cela ne rend pas compte de la cinématique du mouvement. 


e On peut en tout cas faire figurer la « flèche du temps » à l'extré- 
mité de cette trajectoire, qui indique le sens de parcours lorsque le 
temps croît. Cette façon de faire peut être reproduite sur un même 
dessin pour quelques valeurs de X,, suffisamment diverses pour 
être significatives des différentes trajectoires. C'est ce que propose 
la figure 1 ; sur une même orbite, les positions qui sont du même 
côté de X que la flèche du temps sont dans l'avenir de X, les autres 
sont dans son passé. 


Figure 1 - Flèche du temps 
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Figure 2 - Éléments de contact 


e Une autre façon de procéder est d'introduire le champ de vec- 
teurs X AX. On le visualise ordinairement par un vecteur dont 
l’origine est placée à la position X. Pratiquement, on ne place, bien 
entendu, que quelques-uns de ces vecteurs ! Puisque X” (t) dirige la 
tangente à la trajectoire t- X(t) (ce n'est rien d'autre que sa 
vitesse), l'égalité X’(t) = AX(t) exprime que, en chaque point 
d'une trajectoire, la vitesse est égale au vecteur du champ en ce 
point. Compte tenu que la valeur exacte de X” (ft) nous intéresse 
moins que sa direction, et aussi pour des questions de lisibilité, 
nous ne visualiserons pas le champ lui-même, mais seulement, par 
des fléchettes, la direction de AX. On parle alors d'éléments de con- 
tact. Bien entendu, lorsque AX = 0, on ne fait rien figurer (la direc- 
tion n’est pas définie). 


e La figure 2 propose à la fois des éléments de contacts et des 
orbites. Nous utiliserons cette présentation dans la suite. On dit que 
l'on a visualisé le portrait de phase de l'équation différentielle 
X' = AX. Cette possibilité est offerte pour de petites valeurs de p: 


— comme ici, lorsque p = 2 ; 
— éventuellement, par un dessin perspectif, lorsque p = 3. 


Le cas, très simple, où p=1, sera étudié en détail dans le 
paragraphe 5.1. 


L'espace RP dans lequel arrivent les solutions t- @(t, X;) de 
l'équation différentielle est appelé espace des phases. 


1.2 Exemple d'une équation scalaire 
d'ordre deux, linéaire 
et à coefficients constants 


Considérons l'équation différentielle scalaire, d'ordre 2, linéaire 
et à coefficients constants : 


(F) | x"+x = 0. 
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Ses solutions sont, par définition, les fonctions x, deux fois dériva- 
bles, de R, vers R,, satisfaisant à l'égalité : 


VteR x”(t)+x(t) = 0. 


Bien que nous puissions parfaitement résoudre cette équation, 
nous ne le ferons pas dans un premier temps. 


Si xest une solution, introduisons : 
X(E) = (x(P, x (©). 
Nous voyons que, alors, X vérifie : 


X'(E) = (x, x”) = (XP, -x(P)), 


0 1 


donc, en posant À = | ] et en identifiant les vecteurs de R? 
aux matrices colonnes qui les repèrent dans la base canonique, 
que : 


X'(t) = AX(t). 


Réciproquement, si X= (x, y) est solution de l'équation différen- 
tielle vectorielle de taille 2 et d'ordre 1: 


(E) X' = AX, 


on a: 
y=x et x”+x= 0. 


Il y a donc une bijection entre l'ensemble des solutions x de 
l'équation scalaire (F) d'ordre 2 initiale et l'ensemble des solutions X 
de l'équation vectorielle (E) d'ordre 1, mais de taille 2. Cette bijection 
est l'application : 


x (x, X°). 


Bien que, dans l'optique d’une étude pratique, il soit souvent con- 
cevable de conserver l'équation scalaire x”+ x = O0 sous sa forme 
initiale, le point de vue théorique réclame la mise sous la forme 
équivalente du système d'ordre 1 (mais de taille 2) X’ = AX. Pour 
le comprendre, il faut mesurer ce que doit être l’espace des phases. 
Si l’on conserve l'équation du second ordre x”+x = 0, il seraït 
naturel de considérer que l'espace des phases est l’espace d'arrivée 
de l'application x, c'est-à-dire R . Pourtant, une solution est détermi- 
née par le couple (x{0), x'(0)}, et non pas seulement par la valeur 
x(0). Cette situation correspond à beaucoup de configurations de la 
mécanique, où la loi du mouvement est de la forme x” = fx x’), 
auquel cas l'avenir et le passé du point matériel sont déterminés, 
non seulement par sa position « géographique » à l'instant initial, 
mais encore par sa vitesse à cet instant : autrement dit, par la valeur 
(x (0), x’ (0)), qui correspond bien à un point de R2. Nous revien- 
drons amplement sur de tels exemples dans la suite. 


Dans ces conditions, l'espace des phases est l’espace d'arrivée de 
X, soit R?, que l’on désignera aussi par l'espace des (x, x") pour en 
rappeler l'origine. La nomenclature introduite dans le para- 
graphe 1.1 s'applique sans changement. Pour tracer le portrait de 
phase, nous pouvons, dans le cas qui nous occupe, résoudre expli- 
citement l'équation différentielle. On obtient alors l'expression de 
X: 


X(E) = (x(0)cost+ x’(0)sint,-x(0)sint+x’(0)cost), 


soit, en termes de flot : 
CE Xo) = (Xycost+ x,sint, -x;sint+ x, cost) 


lorsque X5 = (X3, X2). 
Les orbites sont alors faciles à décrire. Si l'on munit R2? de sa 
structure euclidienne canonique, il s’agit : 


— de cercles centrés en (0,0), lorsque (x;, x,)# (0, 0) ; 
— de l'origine lorsque (x, x,) = (0, 0). 
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Figure 3 - Portrait de phase dans le plan des (x, x) 


La figure 3 rend compte du portrait de phase. 


Il faut bien noter que le champ d'éléments de contacts donne 
la direction de la « vitesse » dans l’espace des phases, c'est-à- 
dire la direction de X”, mais n’a rien à voir avec la direction de 


la vitesse du mouvement tr x(t) (qui d’ailleurs, dans le cas 
présent, est à valeurs dans la droite réelle). En fait, notre étude 
se situera, pour l'essentiel, dans l’espace des phases, et non pas 
dans l’espace « géographique ». 


1.3 Exemple d'une équation scalaire 
d'ordre deux, non linéaire, autonome 


Considérons l'équation différentielle : 
(F) x"+x + sinx = 0. 


On dit qu’elle est autonome ; cela signifie qu’elle est de la forme 
x” = f(x X), par opposition à l'équation différentielle scalaire 
d'ordre 2 la plus générale x” = g(t, x, x’), où la loi qui régit l'évolu- 
tion du phénomène dépend du temps. Bien entendu, le fait qu'une 
équation différentielle autonome traduise une loi indépendante de t 
ne signifie nullement que ses solutions, elles, soient constantes au 
cours du temps. L'exemple du paragraphe 1.2, dans lequel l'équa- 
tion est autonome (et en outre linéaire), doit nous en convaincre. 


Cette équation différentielle autonome d'ordre 2 équivaut à un 
système vectoriel autonome de taille 2, et d'ordre 1. La démarche 
est analogue à celle suivie dans le paragraphe 1.2. En posant 
X = (x x), on voit que X'est solution de : 


X' = F(X) avec F(x3, x2) = (2%, -sinx3 - x) 


et réciproquement. 


L'équation différentielle vectorielle (ou encore système différen- 
tiel) X’ = F(X) est naturellement qualifiée elle aussi d'autonome. 


Dans le cas présent, il n’est pas possible d’expliciter le flot à l'aide 
des seules fonctions élémentaires, ni même en y adjoignant le sym- 
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Figure 4 - Portrait de phase du système (F) 


bole / de primitivation. Cela ne signifie pas que l'on ne sait pas étu- 
dier le flot. Le portrait de phase de la figure 4 est obtenu à l'aide de 
méthodes numériques de résolution approchée. 


2. Outils mathématiques 


Dans ce paragraphe 2, nous introduisons le matériel relatif aux 
équations différentielles. Nous donnons de manière formelle le lan- 
gage des systèmes dynamiques dans le cadre du flot d'une équation 
différentielle autonome. 


2.1 Théorème de Cauchy et Lipschitz 


Définition 1. 

Soient £ un R-espace vectoriel de dimension finie p, Q un 
ouvert non vide de Rx E et fune application de Q vers E. On 
appelle solution de l'équation différentielle d'ordre 1 et de taille 


p: 
(F) KE MEN) 


toute application X, de source un intervalle Z de R, de but E, 
dérivable sur I, telle que, pour tout te 1, (t, X(t)) appartienne à 
Q, et telle que, de plus : 


Vtel X'(t = f{t X(t). 


Une telle solution sera désignée par (1, X). 
Soit, en outre, (t5, X)) un élément de Q, appelé condition ini- 
tiale. On appelle solution au problème de Cauchy :_ 
X°(E) = FE X(E)) ; 


toute solution (1, X) de l'équation (F) telle que, de plus, t appar- 
tienne à J et telle que X(t) = X5. 


Le langage suivant est d'usage. La variable't est appelée la varia- 
ble temporelle ; ty est appelé instant initial, X; valeur initiale, 


1 intervalle de temps. Bien entendu, ce langage peut varier en fonc- 
tion des situations modélisées. Nous verrons, dans la suite, que 
nous emprunterons en général la terminologie à la physique et, plus 
précisément, à la mécanique classique. 


Notons S((F), to, Xo) l'ensemble des solutions du problème de 
Cauchy précédent. | 


Si (1, X) appartient à S ((F), t, Xo), il est facile d'obtenir d’autres 
solutions : en restreignant X à un sous-intervalle J de 7, avec la con- 
trainte que J contienne t. Il est clair que la restriction Y de X à J défi- 
nit encore un élément (J, Y) de S{((F), to, Xo). On dira, dans ces 
conditions, que (X, 1) prolonge (Y J), et on notera cette situation : 


(Y, J) < (XD) 


Le symbole <, qui permet de comparer deux applications, ne doit 
pas être confondu avec le symbole <, qui représente l'ordre usuel 
dans R. 


Il est très facile de constater que < est une relation d'ordre, non 
total. Considérons le plus grand élément (Z, X) de S{((F), to, Xo) 
(pour la relation < ), s'il existe. Il vérifie la propriété suivante. 


Tout élément (J,Y) de S ((F), to, Xo) est la restriction de (I, X), ce 
qui signifie que JCI et que, pour tout t dans J': 


Y( = X(D. 


Théorème 1 (de Cauchy et Lipschitz). 


Avec les notations de la définition 1, supposons f de classe €! 
sur Q. Il existe alors un unique plus grand élément (1, X) de 
S((F), to Xo). De plus, l'intervalle 7 est ouvert. 


En outre, soit € un compact de Q. Pour tout f, € 1 : 
— il existe t =t, telque (4 X(f))e C ; 
— de même, ilexiste t <t, tel que (#, X(f))e C. 


En particulier, (X,1) est certainement un élément maximal de 
l'ensemble des solutions de l'équation (F), ce qui signifie que l’on ne 
peut prolonger X en une solution de l'équation (F) à un intervalle 
strictement plus grand que 7. Réciproquement, un élément maxi- 
mal (4, X) de l’ensemble des solutions de l'équation (F) est plus 
grand élément de S((F), ti, X (t)) pour tout #, de 7. De ce fait, on 
appelle souvent une telle solution une solution maximale (bien 
qu'elle vérifie une propriété meilleure). 


# Le théorème de Cauchy et Lipschitz rend compte d’une observa- 
tion rendue naturelle par de nombreux exemples issus des sciences 
de la nature. Lorsque la loi d'évolution temporelle d’un système -— loi 
résumée ici par l'application f, qui fournit l'équation différentielle 
régissant le phénomène -— est connue, et est d'ordre 1, et lorsque le 
système est entièrement déterminé à un instant donné t (par la 
valeur Xo), l'avenir et le passé du système sont entièrement détermi- 
nés. En vérité, cette constatation doit être approfondie à la lumière 
du théorème 1. 


Tout d’abord, il est bien nécessaire que l’évolution du système 
soit régie par une équation différentielle. D'autre part, cette équation 
différentielle ne peut être absolument quelconque : si fest simple- 
ment supposée continue, le théorème est en défaut. Cette remarque 
met en évidence qu'une loi physique « déterministe » exprimée par 
une équation différentielle doit satisfaire certaines conditions de 
régularité. 


Un autre aspect du théorème est la nature de J7, qui reste en 
général inconnu. Il n'est en particulier nullement obligatoire que 7 
soit égal à R, même lorsque Q = RXE. Nous indiquerons la 
dépendance de 1 par rapport aux conditions initiales en le notant, si 
besoin est, Z(t,, X,ÿ). Nous savons néanmoins que 7 est ouvert. 
Pour le déterminer, nous utiliserons souvent la dernière partie du 
théorème, qui peut s'exprimer en disant que le graphe de (I, X), 
dans tout avenir, quitte tout compact de Q. Cela ne nous dit rien lors- 
que sup] = +. En revanche, supposons que sup] = b soit un 
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réel. Nous pourrons exploiter le résultat comme dans l'exemple sui- 
vant. 


Exemple : supposons que Q = RXE et qu'une solution X soit 
bornée sur 1. Autrement dit, X(t) appartient à une certaine boule B de 
E pour tout t. 

Montrons qu'alors Z = R. Dans le cas contraire, Z admet une borne 
réelle ; supposons par exemple que c'est la borne supérieure b. En par- 

. ticulier, il existe un sous-ensemble compact [t,, blx B de RxE con- 
tenant le graphe de X dans l'avenir de t,, ce qui est une contradiction. 


Le théorème a une conséquence immédiate. Soient (J;, Y,) et 
(J, Y:) deux éléments de S ((F), to, X0). Ils sont prolongeables par 
(I, X). En particulier, si : 


te Jin, Yi = Y( = XP. 


On exprime ceci par l’unicité locale des solutions du problème de 
Cauchy. 


Désormais, nous ne nous intéresserons qu’à la plus grande solu- 
tion du problème de Cauchy. Puisqu'elle dépend de la condition ini- 
tiale, nous poserons : 


X(t) = 4 (6 ty Xo): 


2.2 Régularité par rapport 
aux conditions initiales 


Nous nous plaçons, avec les notations du paragraphe 2.1, sous les 
hypothèses du théorème 1. Nous disposons, pour chaque point 
(to Xo) de Q, de la plus grande solution du problème de Cauchy : 


It X0) — E 
tr CE Lo X0) 


Dans ces conditions, notons 4 l'ensemble des (# t5, X,) tels que 
te Ito Xo). L'application $ est définie sur A. Supposons f de 


classe CK, kétant dans [1, ++]. Il est clair que : 
te O(6 ty Xo) = X(f) 


est de classe Ck+1, 


En effet, X est dérivable sur Z, donc continue. Comme 
X'(t) = f(t, X(H)), X' est continue, donc X est de classe C?, La 
même relation prouve alors que X est de classe C?, et une récur- 
rence immédiate termine le raisonnement. 


Cependant, nous ne connaissons rien a priori de la dépendance 
de par rapport aux conditions initiales. Le théorème 2 répond à 
cette question. 


Théorème 2. 


Avec les notations qui précèdent, supposons f de classe CE &k 
étant dans [1, + ].L'ensemble % est un ouvert de RXRXE. 
De plus, l'application & est de classe CKsur 4. 


Ce théorème exprime essentiellement que, si l'on se donne deux 
conditions initiales proches, et deux instants proches en lesquels les 
solutions sont respectivement définies, les valeurs des solutions à 
un instant donné sont proches (continuité de #); de façon plus 
générale, la dépendance de la solution par rapport à l'instant et aux 
conditions initiales est de même classe que f. 
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2.3 Cas des équations non scalaires 


Une équation différentielle très générale, de la forme : 
(G) x) = DER os MIN, 


où gest une application d'un ouvert Q de R x E9, est dite d'ordre q,. 
Ses solutions sont les applications x, q fois dérivables sur un 
intervalle 1 de R,, à valeurs dans E, telles que, pour tout t dans 1, 
(x x, …, xX4-1)) appartienne à Q, et telles que, de plus : 


Vter x0)(t) = g(t,x(8, x'(t), …, xt4- D (E)). 


On se ramène au cas d’une équation scalaire d'ordre 1 en consi- 
dérant l'application X, dérivable sur 7 et à valeurs dans E', définie 
par : 


XP = CE, XP, …., x DD) ; 


si xest solution de l'équation (G), X est alors solution de l'équation 
d'ordre 1 : 


(F) X' =f({tX), 


où : 
K(t, Xo; …. XQ-1) = (X3; .…. Xq-1 g(t, XO» es 3 Xq-1)) . 


Ainsi, fest définie sur Q et est de même classe que g. Réciproque- 
ment, si X est solution de l'équation (F), sa première composante 
sera solution de l'équation (G). Cette bijection entre ensemble des : 
solutions de l'équation (F) et ensemble des solutions de l'équation 
(G) induit une bijection entre S ((F), to X0) et une partie de l’ensem- 
ble des solutions de l'équation (G) ; précisément, si 


X0 = CX0, o& Xo, 1 +++ X0, g-1) 7 
la condition X(t,) = X, équivaut à : 
X(t0) = XG, 0» ….) x(a- (4) = X0, q-1 « 


Ainsi, nous devons, pour une équation d'ordre q, formuler le pro- 
blème de Cauchy par un q-uplet de conditions initiales : 


X(to) = X0, 0 ne Ce = X0, q-1 : 


Nous noterons encore S ((G), t, Xo) l'ensemble des solutions de 
l'équation (G) vérifiant cette condition. 


Nous avons pu ramener les équations d'ordre q aux équations 
d'ordre 1. La relation < sur les applications étant naturellement 
analogue à celle introduite au paragraphe 2.1, les théorèmes 1 et 2 
se paraphrasent sans difficulté, comme dans le théorème 3. 


Théorème 3 (de Cauchy et Lipschitz). 

Avec les notations précédentes, supposons g de classe C\ sur 
Q. Il existe alors un unique plus grand élément (ZX) de 
S ((équation (G)), t, Xo). De plus, l'intervalle J est ouvert. 

En outre, soit C un compact de Q. Pour tout f; € I : 

— il existe tzt, telque X(téE C ; 

— de même, il existe t < t, tel que X(F)# C. 


Cet énoncé est conforme à la situation physique la plus classique ; 
lorsque l'accélération est fonction de la position et de la vitesse (et 
éventuellement du temps), et sous réserve que cette dépendance 
soit assez régulière, la connaissance, à un instant t,, de la position 
x (to) et de la vitesse x” (t,) permet de déterminer l’évolution passée 
et future du système. 
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2.4 Équations autonomes 


Définition 2. 

Soient @ un ouvert de E et f une application de Q vers E. 
L'équation différentielle : 
(E) X'=f( 


est dite autonome d'ordre 1. 


Nous sommes dans la situation générale évoquée dans le 
paragraphe 2.1, disposant en outre de l'information supplémentaire 
que f dépend seulement de X (et non pas du couple (t, X)). Comme 
on l’a vu dans le paragraphe 1.1, on peut interpréter f comme un 
champ de vecteurs sur Q. L'espace E est souvent appelé espace des 
phases. En fait, f n'est définie que sur une partie de E. Il est donc 
plus légitime de considérer l'ouvert Q, que l'on appelle la variété des 
phases. 


Nous noterons qu’une équation d'ordre q, de la forme : 
xD = g(x x, .…., xD), 


que l'on appelle équation différentielle autonome d'ordre q, se 
ramène à la forme précédente par la méthode du paragraphe 1.3, ce 
qui justifie que, dans la suite, nous étudiions exclusivement l'équa- 
tion (E). 


L'ensemble des résultats énoncés dans le cas général s'applique, 
avec quelques précisions. Nous nous placerons désormais dans le 
cas où fest de classe C!, ce qui permet d'appliquer, en particulier, le 
théorème 1 de Cauchy et Lipschitz. Tout d’abord, nous noterons la 
propriété d’invariance par translation temporelle suivante. 


Proposition 1. 
Soient t, et À deux réels. Pour (I, X)e S ((E), to X0) , PoOsons : 


Vuel+2 Xj(u) = X(u-À). 
L'application : 
S ((E), ty X0) > S (E) to + À, Xo 
XF X; 


est une bijection. 
Preuve. © Soit (1, X)e S ((E), to Xo)-: On a: 


Vtel X'(t) = AX(#)). 


Soit ue 1+À, où 1+À désigne l'intervalle obtenu à partir de 7 
par translation de À. Puisque u = t+1 avec te 1,0na: 


Xi Qu) = X'(u-A) = X(8 = FX) = FOQU-A)) = FO CU)). 


Donc X; est solution de l'équation (E). 


De plus : 
Xito + À) = X(to) = Xo: 
Donc : 
Xe S ((E), to + À, Xo) - 
Le caractère bijectif de l'application est alors immédiat. © 


La proposition 1 a pour conséquence immédiate que les plus 
grandes solutions (au sens de la relation < ) se correspondent elles 
aussi par translation. Précisément, on a: 


CE ti Xo) = O(t+(t1- to), os Xo) 
et It Xo) = (to Xo) + (t1— to) 
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Cette invariance par translation temporelle permet, sans restric- 
tion, de fixer arbitrairement un instant initial, que l’on prend, dans 
toute la suite, égal à 0. En d'autres termes, on ne s'intéressera 
qu'aux solutions définies en l'instant 0, les autres pouvant être obte- 
nues, si on le désire, par le moyen précédent. 

Nous pouvons alors simplifier les notations. On notera I(X,) au 
lieu de (0, X,) et, de même, ÿ(& Xj) au lieu de ÿ(t, 0, X,). Ainsi, 
% est l'ensemble des couples (t, Xj)e RXQ tels que tappartienne 
à 1(X) . D'après le théorème 2, X est un ouvert de RXxE.En outre, 
si le champ de vecteurs f est de classe CF, l'application est de 
classe CK 


Définition 3. 
L'application #, définie sur l’ouvert , est appelée flot du 
champ de vecteurs f (ou de l’équation différentielle X'=f(X)). 


Pour comprendre cette définition, on peut considérer le champ de 
vecteurs comme un champ d’impulsions ; une particule occupant à 
l'instant 0 la position X, (dans l'espace des phases) occupera à l'ins- 
tant t la position #(t, X,). En d'autres termes, soumises à ce champ 
d'impulsions, les particules de l'espace des phases s'écoulent selon 
le flot à. 


Posons, pour (t, X)e A, p(Xo) = O(E Xo)- Ainsi, t-> p(Xo) est 
la solution de l'équation (E) qui, à l'instant initial, prend la valeur X0- 

Proposition 2. 

On a l'égalité : 

(Pr PSN Xo) = Pt sX0) 

lorsque, d'une part p%XQ) et, d'autre part, @,4(Xo) ou 
(#, o p)(X) sont bien définis. 

De plus : 

Preuve. © La deuxième égalité résulte de la définition, puisque 


#(0, X,) est la valeur à l'instant initial de la solution prenant à l'ins- 
tant initial la valeur X6. 


D'un autre côté, tr p,, %XQ) est solution maximale de l'équa- 
tion (E) sur 1(X,)-s grâce à l'invariance par translation temporelle, 
et prend la valeur p,(X,) à l'instant 0 si 0e I(X)-5, c'est-à-dire si 
se I(Xo)- 


Or tr p{pXo)) est la solution maximale de l'équation (E) qui, 
à l'instant initial, prend la valeur p.(X5) . 


Ces deux applications ont donc même intervalle de définition, et y 
coïncident. Ô© 


Corollaire. 
Soient te R et A, l’ensemble des X, tels que (t, X0) € À. L'appli- 


cation y, est un difféomorphisme de classe C 1 de U, sur U,. 


Preuve. © Puisque 4 est ouvert, U, est aussi ouvert. De plus, si 
Xe Ur, P(Xo) est bien défini, et p;_4(Xo) = X, l'est aussi. Donc : 


Ida, 


P_+0O + 
De même : 
Pr0 P+= Ida, 


Comme et +, sont de classe CK+1, le résultat en découle. © 
t À Pr 


Définition 4. 

Les arcs paramétrés tr g(t Xi), pour XQ dans Q, sont appe- 
lés trajectoires du système dynamique. L'image d'une telle tra- 
jectoire est appelée orbite du système dynamique. 
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Soit XQ un point de Q. Il existe, d’après le théorème de Cauchy et 
Lipschitz, une orbite contenant le point X4, à savoir l'image de 


t=@(t X5) = X(f). Supposons que X; appartienne à l'image de 
tRO(E X;) = Y(). Cela signifie qu'il existe t, dans Z(X,) tel que 
X(0) = Y{f,). La proposition 1 entraîne alors que Y(t) = X(t-t:;) 
pour tout t dans Z(X;). 


Une première conséquence de cela est que l’image de Y'est égale 
à celle de X, donc que X, appartient à une seule orbite, que l'on 
appellera l'orbite de X5. En particulier, les orbites réalisent une par- 
tition de Q. 


En second lieu, l'application X peut être injective sur 7, ce qui 
signifie que l'orbite n’est parcourue qu’une seule fois, ou encore 
que, pour tout point M de l'image de X, il existe un instant t, et un 
seul, tel que X(t) = M. 


Dans le cas où X n’est pas injective, il existe deux instants t; <f, 
tels que : 


X(t1) — X(b), 
et donc, pour tout t dans I1,ona: 
X() = X(t+ 2) 


avec À = tete 
En d'autres termes, X est une application périodique. 


Si X n'est pas une application constante, soit 7 sa plus petite 
période strictement positive. Ce qui précède montre que, sur [0, TT, 
X est injective, et que X(0) = X(7). La trajectoire est dite fermée 
simple. On notera que l'adjectif « fermé » n'a rien à voir avec celui 
utilisé en topologie. Dans le contexte présent, il indique simplement 
que la courbe se referme sur elle-même. 


Si X est constante (égale à Xi), on a: 
X(D =0, 


donc f(X) = 0. 


Réciproquement, si X, est un élément de Q tel que f(X,) = 0, 
l'application constante égale à X, sur R est solution de l'équation 
différentielle. Une telle orbite, réduite à un point, est appelée posi- 
tion d'équilibre. 

Soit X une solution non constante. Puisque les orbites sont dis- 
jointes, son image ne peut passer par une position d'équilibre. En 
particulier, puisque X’(t) = f(X(t)), on ne peut avoir X’(t) = 0. 
Ainsi, la trajectoire X est un arc régulier, c'est-à-dire qui admet une 
tangente en tout point. 


Résumons. 
Proposition 3. 


Soit X° = f(X) une équation autonome, où f est de classe CK sur 
Q. Les orbites réalisent une partition de Q. Une trajectoire est : 


— ou bien une application constante sur R , qui correspond à une 
orbite réduite à un point, appelée position d'équilibre ; 

— ou bien une application périodique sur R , qui correspond à un 
arc de classe C **T fermé, régulier, simple, d'image compacte ; 

— ou bien une application injective sur l'intervalle ouvert I, qui 
correspond à un arc de classe C F5} régulier, simple, d'image non 
compacte. 


Preuve. © Il ne nous reste qu'à prouver la non-compacité d’une 
trajectoire injective. 


Supposons par l'absurde que X, solution de l'équation (E), est 
injective et que X(1) est compact. Montrons d’abord que 7 = R :si 
l’on suppose par exemple b = supl< + ,{(t, X(f)) reste dans 
[0, bJx X(1) pour t> 0, ce qui contredit le théorème de Cauchy et 
Lipschitz. 

Comme image d’un compact par l'application continue X, X([-n, nl) 
est un compact, donc un fermé. 


Montrons que son complémentaire est dense dans X(R). Pour 
cela, introduisons l'ensemble À des limites des suites (X(f,)), lors- 
que (f,) est une suite quelconque de réels telle que (éd) tende vers 


+, Puisque X(R) est compact, A est inclus dans X(R). En outre, 
A est non vide : il suffit de considérer n'importe quelle suite de réels 
(|ud) tendant vers +, et en extraire, par compacité de X(R) , une 


suite (|v4) telle que (X(v,)) tende vers a. Mais, quel que soit f : 
Pt (X(0)) = PP (X(O))) . 


Puisque Py(X(0)) = X(vV,) = &, la continuité de 9, entraîne 
que : : 


Pta (X(0)) — pe). 


Par conséquent, p,(&) appartient à À. Il en résulte que A contient 
tous les éléments de l'orbite de æ& et donc contient X(R). 


Soit alors fun élément de X(R) ; il est dans À. En particulier, pour 
tout n, il existe une suite (|t) de réels tous de valeur absolue plus 
grande que n, telle que X(t,) — B. Cela entraîne que X(R —-[-n, n]) 
est dense dans X(R). Comme X est injective, X(R —-[-n, n]) n'est 
autre que le complémentaire, dans X(R), de X([-n, n]). Par consé- 


quent, X(R-[-n, n]) est un ouvert dense dans le compact X(R).Le 
théorème de Baire entraîne alors que l'intersection de tous ces 
ouverts (lorsque n décrit N ) contient au moins un point, de la forme 


X(t), mais aussi de la forme X(t) pour un t’ tel que |t{>n>l{t|. 
Cela contredit l’injectivité de X. © 


2.5 Retour sur le cas non autonome 


Considérons à nouveau une équation non autonome : 
(F) x =f(t x). 


Posons X(t) = (t, x(t)). Il est facile de constater que X est solu- 
tion de l'équation autonome : 


(E) X° = g(X) 


avec g(Xx4, X2) = (1, F(X7, X2)). 

Cependant, la réciproque n'est pas vraie, puisque, si 
X(#) = (x(P,xX(t)) satisfait l'équation (F), la primitivation de 
x, = 1 fournit seulement : 


XA(u) = u+c. 
Si on se penche sur le problème de Cauchy : 
ATLAS » A0 = Xe: 
on obtient la condition supplémentaire : 
X(0) = (0, x(0)), c’est-à-dire x,(0) = O0 et x,(0) = x, 


qui rend cette fois équivalents les deux problèmes. Ainsi, étudier 
une équation non autonome avec une condition initiale (0, xo) 
revient à étudier le flot d’une équation autonome limité à des 
Valeurs de X, de première composante nulle. Malheureusement, 
dès que la condition initiale est du type x{(f,) = x, avec t, non nul, 
on ne peut en général se contenter d'étudier le flot du système auto- 
nome précédent, puisqu'il n’est nullement certain que (#5, xo) soit 
dans 4. On gardera donc à l'esprit qu'une équation non autonome 
peut se ramener à une équation autonome, assortie des conditions 
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initiales ad hoc, mais que la théorie des systèmes dynamiques, et 
notamment l'utilisation du flot, ne s'applique pas sans restriction. 


Exemple : soit (F) l'équation scalaire x'= 2X2t, choisie parce que 
l'on peut en expliciter les solutions. Elle est à variables séparables. Sa 
résolution conduit aux différentes solutions suivantes : 


I=RxÔ0=0; 


1 


I1=R,xX( 0 = —-———, 
es z 


a>0 ; 
21 En es 0 : 
I = —c, —a[, X( ) . ne ' 


1 , 
1 = ]-a, al, xX(Ô = n'es ; 


4 
1 =1]a, +œl[,x(t) = part. 


Les trajectoires correspondant à ces équations non autonomes 


sont tracées sur la figure 5. 


D'un autre côté, si l’on étudie le flot, on constate que, dans 
l'espace des phases, on n'obtient comme solutions définies en 0 que 
celles dont le graphe est tracé en gras (non tireté). 


3. L'exemple du pendule 


Nous allons dans ce paragraphe introduire un système, approxi- 
mant un système physique élémentaire. Après avoir modélisé ce 
système, nous constaterons que son évolution obéit à une équation 
différentielle autonome d'ordre 2, c'est-à-dire à un système différen- 
tiel autonome de taille 2. Nous étudierons alors les solutions, dans 
le but à la fois d'utiliser les outils introduits dans le paragraphe 2, de 
tracer le portrait de phase et d'illustrer le langage propre aux systè- 
mes autonomes. Nous démontrerons aussi, volontairement, des 
résultats qui, dans ce cas, pourraient paraître évidents (qui n'a pas 
assisté à cette expérience, au moins virtuellement ?), principale- 
ment pour pouvoir les appliquer plus tard dans des contextes où le 
but de la théorie serait de prévoir, et non pas de confirmer, des résul- 
tats physiques. 


3.1 Le problème physique 
et sa loi d'évolution 


On considère une tige rigide de longueur £ , à l'extrémité E de 


laquelle est fixée une masse m. L'autre extrémité du pendule peut 
tourner autour d’un point fixe F ; la masse de la tige est supposée 
nulle. Le pendule est placé dans un champ gravitationnel constant 
de l’espace affine euclidien de dimension 3. Constant signifie ici 
indépendant de la position et de l'instant. On supposera, pour con- 
server une situation traditionnelle, que le champ est d'intensité g, de 


> > 
sorte que l'accélération subie par une masse mest mg i , où i est 
unitaire. On compiète ce vecteur en une base orthonormale 


> > > 
(i{,j,k) de l'espace. Nous supposerons en outre (c'est une restric- 
tion importante au cas général) que le pendule se déplace dans le 


> > 

plan de repère (F,i,j) plan vertical contenant F. Le point F se 
déplace alors sur un cercle. Dans ce plan orienté, on introduit un 
angle polaire 8, comme sur la figure 6. : 


A ce sujet, il faut constater que la position de.E n’est pas bijective- 
ment liée à la valeur de 6, définie à un multiple entier de 2x près. Ce 


Figure 5 - Flot d'une équation non autonome 


Figure 6 - Pendule simple 


faisant, nous avons fait une hypothèse implicite sur la variété des 
phases, en paramétrant le cercle par 6. Nous n'étudierons pas ici la 
validité de cette hypothèse. Disons simplement que, de cette façon, 
nous pourrons nous placer dans le cadre de cet article. 


L'accélération tangentielle à laquelle est soumis le point E est 
égale à: 


On obtient donc : 
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Cette équation différentielle autonome et scalaire d'ordre 2 régit 
l’évolution du système. 


Posons : 


On obtient : 


w "(#) = 207) 5 -sin(d) = _sin(y(#). 


Ainsi, nous avons pu nous ramener au cas où a= 1, par un chan- 
gement de variable temporelle qui, ici, correspond à prendre sim- 
plement une autre unité de temps. Ce procédé, de normalisation, est 
couramment utilisé pour alléger la présence des paramètres. En 
revanche, lorsque nous voudrons étudier des systèmes dynamiques 
dépendant d'un paramètre, nous ne pourrons recourir à cet expé- 
dient qu'avec parcimonie, puisque le changement de variable 
dépend justement du paramètre. 


Nous étudions à présent l'équation différentielle : 


(F) 0"= -sin6. 


3.2 Utilisation d'une intégrale première 


L'équation (F) équivaut au système : 
(8, 8) = f(6, 0”) 


avec f(x, X2) = (X2,-sinx;). 

L'espace des phases est, ici, R?, que l’on appellera, pour d'évi- 
dentes raisons, l’espace des (0, 0’). A l'occasion, on désigne même 
par (8, 8‘) un élément de R2. 


3.2.1 Premières remarques 


B Étude de 7 


Soit 8 une solution maximale, sur 7, de l'équation (F). Supposons, 
par l'absurde, que Z admet la borne supérieure réelle a. Puisque 
9” = -sin0, 0” est bornée sur 7, donc a une intégrale convergente 
en b-. Cela implique que 0” admet une limite en b7 puis, par un 
argument analogue, que 6 admet une limite en b-. Cela entraîne 
que 8 est bornée sur [0, b[, donc que le graphe de 8 reste dans un 
compact pour {> 0, ce qui est une contradiction. 


Ainsi, Z n'est ni majoré ni, par un argument identique, minoré. 
Donc : 
I=R. 


E Recherche des positions d'équilibre. 
Les positions d'équilibre sont les solutions de l'équation : 


(x, -sinx,) = (0,0). 


On obtient les positions (kx, 0) (pour ke Z ). Physiquement, la tige 
du pendule est alors en position verticale (vers le bas, si k est pair, 
vers le haut sinon), et le point E a une vitesse nulle. 


Dans le paragraphe 3.2.2, on s'intéresse aux trajectoires non 
constantes. 


3.2.2 Obtention d'une intégrale première 


Multiplions par 24’ l'égalité : 
0”+sin0 = 0 
et primitivons l'égalité ainsi obtenue. || vient : 
(G) 9’2-2cos6 = c. 


On dit que l’on a obtenu une intégrale première du système. lci, 
l'intégrale première est l'application £ définie sur R2 par: 


E(X1, X2) = X3-20c0Sx1. 


De façon générale, on appelle intégrale première d'une équation 
différentielle autonome toute application, de classe C\ sur Q, cons- 
tante sur chaque orbite. 


Des intégrales premières sont dites indépendantes lorsque, pour 
tout X, dans Q, les applications linéaires tangentes en XQ sont indé- 
pendantes (ou de manière équivalente, dans le cas d'une structure 
euclidienne, lorsque les gradients sont indépendants). Dans le cas 
d'une seule intégrale première, cela signifie que l'application 
linéaire tangente n'est jamais l'application nulle, ce qui exclut en 
particulier les applications constantes. 


Couramment, « famille d'intégrales premières » est souvent uti- 
lisé pour un raccourci de «famille d'intégrales premières 
indépendantes ». Par exemple, puisque : 


grade = (2sinx, 2X), 


le gradient ne s’annule qu’en les positions d'équilibre ; £ est en fait 
une intégrale première « indépendante » sur R? privé des positions 
d'équilibre. 

Nous savons qu’une orbite est incluse dans un ensemble : 


E(X7, X2) = C, 


courbe de niveau de l'application €. 
Soit réciproquement : 


te (a), x1 (D) 


un arc de classe C!, régulier, dont l’image est incluse dans un tel 
ensemble. Par dérivation, on obtient : 


xi (HCX (+ sinx(#) = 0. 


Les points ten lesquels x; s’annule forment un ensemble d'inté- 
rieur vide, faute de quoi x, serait constante sur un intervalle, et l'arc 
ne serait pas régulier. L'égalité : 


x(t) + sinx,(t) = 0, 


vraie sur un ensemble dense, est donc vraie partout, et x, est solu- 
tion (peut-être pas maximale) de l'équation différentielle. 


Pratiquement, cette intégrale première va être utilisée pour tracer 
les ensembles : 


E(X7; Xo) = C. 


Il est souvent difficile d'exploiter ces ensembles de façon rigou- 
reuse, puisqu'il faut savoir s'ils peuvent être paramétrés comme ci- 
dessus. Néanmoins, ils donnent des indications précieuses sur le 
portrait de phase, permettent de le visualiser lorsqu'on en a démon- 
tré les propriétés, et de retrouver des propriétés (mathématiques, 
expérimentales ou intuitives) que l’on en connaît déjà. 
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Figure 7 - Courbes de niveau 


La figure 7 donne les courbes de niveau de l'application €. Les 
positions d'équilibre sont indiquées par un petit rond, les courbes 
fermées sont en trait gras cyan et les courbes non compactes en 
tireté cyan. 


On observe tout d'abord que les courbes fermées pourraient 
correspondre à des orbites fermées, donc à des solutions pério- 
diques. Les courbes en tireté laissent penser qu'elles sont asso- 
ciées à des trajectoires non périodiques, mais qui le deviennent 
lorsque l’on identifie x et x1 + 27. 

Enfin, les deux courbes en gras noir ne peuvent correspondre 
qu’à une infinité d'orbites, puisqu'une orbite ne contient pas de 
position d'équilibre. 


3.3 Étude avec conditions initiales 
particulières 


L'étude précédente (5 3.2.2) suggère que chaque orbite coupe une 
droite x = 2kx. Plutôt que de le montrer maintenant, nous allons 
tout d'abord nous limiter à l'étude des solutions 8 telles que : 
8 (0) =0. 


Nous montrerons plus tard qu'imposer cela n'est pas restrictif. 
Nous poserons aussi : 


9’(0) = v. 


On constate aisément que, si 8 est solution de l'équation (F), 
tr 0(-t) l'est aussi. 


On peut donc se restreindre à v > 0. L'équation (G) devient : 
92= 2cos0+v2-2, 


soit : 
(H) g”2= v2-4sin? 2. 


La position de v par rapport à 2 joue un râle évidemment impor- 
tant, que nous mettrons en évidence. 
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E Cas v>2. 


On constate qu’alors 8’ ne s’annule pas, donc que 9' garde un 
signe constant sur R, nécessairement égal au signe de v, c'est-à- 
dire positif. Notons g la primitive de : 


1 


2 in2U 
V 4sin° > 


qui s’annule en 0, c'est-à-dire : 


ur 


L'équation (H) devient : 
g'(8())9"(0 = 1 
soit g(8(t)) = t+ d, avec d= 0 grâce aux conditions initiales. 


Finalement, si h désigne la réciproque de l'application strictement 
monotone g: 


o(t) = h(t). 


On a pu ainsi expliciter la solution 6, comme application récipro- 
que d’une certaine primitive. A l'aide des fonctions de Jacobi, on 
obtient : 


a(t) = jam 2). 


Posons : 


i du 
| 5 _ rein? U 
vi-4sin 5 


g(u+2x) = g(u)+T, 


2 
T = g(27) = ( 
0 
On a clairement : 


et donc: 
h(t+ TD) = h(t) +27. 


Ainsi, entre les temps tet t+ T, le pendule a effectué une rotation 
complète, et se trouve en ces deux instants avec une position (dans 
le plan) et une vitesse égales. On notera que cette « T-périodicité à 
2x près » correspond à une vraie périodicité sur le cercle. 

Cas v < 2. 
Notons toujours g la primitive de : 


1 


fe _ 6 6 
V 4sin* > 


qui s’annule en 0 ; remarquons que g est impaire. D'autre part, g 
n'est définie que sur l'intervalle [- wo, uol, avec : 


ur 


. V 
uQ = 2Arcsins, 


et h, sa fonction réciproque, est définie sur [- to, fol Où : 


Uo 
Le [ du - 
0 [v2-4sin?> 


Sur [- to, tol, on voit, comme dans le cas précédent, que 8 = h. 


Montrons à présent que 6 est 2t-antipériodique. Puisque 
0’(t,) = 0, l'application : 


tr -0(t+2t), 


AF 180 -11 


PRÉSENTATION DES SYSTÈMES DYNAMIQUES 


qui vérifie encore l'équation différentielle (F), vaut en - t&: 
—6(to) = O(-t) ; 
sa dérivée vaut en — to: 
0 (ty) = 0 = 8’(-to). 
Par conséquent, pour tout f, on a: 
O(t) = -8(t+2t). 


Il résulte de cela que 8 est 4t-périodique. 
B Cas v=2. 


Ce cas est particulier, car on peut expliciter la solution à l’aide des 
fonctions élémentaires. On obtient : 


0 
72 _ 2U 
0” = 4cos 5: 


Dans ces conditions, 8” ne peut s’annuler que si 0 = kn, ce qui cor- 
respond à une position d'équilibre, exclue. Donc 8” garde un signe 


0 
constant, nécessairement positif. De même, cos; garde un signe 
positif. Ainsi, l'équation devient : 


0] 
0" = 2cos3, 


à variables séparables, qui s'intègre sous la forme : 


a(t) = 4Arctanet/2,7 1: 


les conditions initiales ayant été utilisées. On constate aisément que 
6 varie entre -x et x. 


3.4 Étude générale 


Justifions à présent que le choix de conditions initiales particuliè- 
res n'est pas restrictif. Montrons pour cela que toute solution 6 dif- 
férente d’un équilibre prend au moins une fois une valeur de la 
forme 2kn. 


Dans le cas contraire, elle prend ses valeurs dans un certain inter- 
valle ]2kx, 2(k+ 1)xl. Quitte à considérer l'application 8 - 2kn, encore 
solution de l'équation (F), on peut supposer que 8 est à valeurs dans 
]J0,2xl. 

H Premier cas. 0” s’annule deux fois au moins dans R. 

Si 0’2t,) = 9’2(t,) = 0, l’application 8”? atteint sur [t, &] un 
maximum, en lequel elle ne peut être nulle (sinon 8 serait constante 
sur [f1, t], et serait donc un équilibre). A cet instant, 288" s'annule, 
donc 8”, puis sin6, sont nuls. 

Mais si 0=x: 

cos 0=-1. 
Or, 0"? atteignant à cet instant un maximum, cos en atteint un 


aussi d’après l'équation (G). Cela entraîne à nouveau que 6 est un 
équilibre. Il ÿ a contradiction. 


# Deuxième cas. 0” ne s’annule pas dans R.. 


Remarquons déjà que 0” =- sin 8 change de signe et s’annule, 
car une fonction convexe (ou concave) bornée sur R est nécessaire- 
ment constante. En un point d'annulation de sin 6, @ prend la 
valeur x. 


Quitte à changer 6 en —- 8, on peut supposer que 0 croît stricte- 
ment. Étant bornée, 0 admet des limites £ et £”’ en -+ et +.1Ilen 
résulte que 8” admet les limites -sin£ et -sin£”, nécessairement 
nulles (sinon |8| tend vers +). Donc: 


€ = O0 et £’ = 2x. 


Mais 9’? tend aussi vers 0 en —-æ et +, ce qui est une contradic- 
tion manifeste avec le tableau de variations de 6 (tableau 1). 


Tableau 1- Variations de 0 (deuxième cas) 


em | à 0 + 
Ce | os + LT 


Æ Troisième cas. 0” s’annule une fois dans R. 


Quitte à changer 6 en — 8, on peut supposer que 6 décroît, puis 
croît. Puisque, comme dans le second cas, 8 admet des limites en 
co et +0, ces limites sont l’un des nombres 0, x ou 2x. 


Puisque 8 prend la valeur r, les limites de 8 en -c et + sont tou- 
tes deux égales à 2x. Il en résulte que 0 ” = - sin0>0 au voisinage 
de -æ, donc que 8” croît au voisinage de -. Comme 8-0 en 
—c , 0 croît au voisinage de -, ce qui est une contradiction. 


Finalement, nous obtenons le résultat voulu : toute orbite coupe 
une droite 0 = 2kn. En utilisant le fait que, si 0 est solution de l’équa- 
tion (F), 0-2kx est aussi solution de l'équation (F), on a bien pu 
ramener l'étude des orbites aux conditions particulières choisies. 


3.5 Interprétation 


Nous avons montré que le phénomène étudié est défini à tout ins- 
tant (8 3.2.1). Lorsque la tige est en position verticale avec une 
vitesse nulle, elle y reste (8 3.2.1). La fonction € (8 3.2.2), qui est une 
intégrale première, peut être interprétée comme l'énergie du sys- 
tème qui, dans notre cas, est la somme de son énergie cinétique et 
de son énergie potentielle. Nous avons montré qu'il existe toujours 
un moment où la tige passe en position verticale, vers le bas, sauf 
dans le cas où elle reste dans la position verticale vers le haut (8 3.4). 
Ce faisant, nous nous sommes ramenés au cas où l'énergie poten- 
tielle est nulle, et nous avons discuté sur l'énergie cinétique du sys- 
tème à cet instant. Lorsque l'énergie cinétique est suffisamment 
grande, le pendule effectue une infinité de tours; dans un mouve- 
ment périodique à 2x près. Lorsque cette énergie cinétique diminue 
pour atteindre une valeur discriminante (v= 2), le pendule effectue 
un seul tour presque complet (en un temps infini). Lorsqu'elle dimi- 
nue encore, le pendule oscille périodiquement, avec une amplitude 
qui diminue elle aussi (8 3). Nous observons une situation peut-être 
étonnante : la position d'équilibre vers le bas, qui correspond à une 
énergie totale nulle, n'est pas asymptotiquement stable au sens où, 
même lorsque l'énergie du système est faible, la trajectoire du pen- 
dule ne tend pas vers cette position pour de grandes valeurs de t. En 
revanche, pour v=2, la position du pendule tend vers la position 
d'équilibre correspondant à une énergie totale égale à 41 


4.1 Le cas « continu » 


Soit une équation autonome : 


(E) X' = f(X) 
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où fest de classe C* sur un ouvert Q de E. Nous disposons d'une 
application #, le flot, de classe CKsur un certain ouvert 4 de RxQ, 
vérifiant les propriétés suivantes : 


pCE CS, Xo)) = p(E+ 8, Xo) ; 
(0, X5) = Xo 


chaque fois que le premier membre de chacune des égalités est 
défini. 

Nous pouvons légitimement définir de façon générale un système 
dynamique à l’aide de ce cas particulier. 


Définition 5. 

On appelle système dynamique « continu » de classe CX la 
donnée d’un ouvert 4 de RXE et d'une application #, de 
classe C* sur %, et de classe CK+1 par rapport à la première 
variable, vérifiant les propriétés suivantes : 


CE DCS, Xo)) = b(E+S, Xo) : 
#(0, X5) = Xo 


chaque fois que le premier membre de chacune des égalités est 
défini. L'application ÿ est appelée flot du système dynamique. 


L'adjectif « continu » mérite une explication. La notion de système 
dynamique que nous venons d'introduire a trait à un flot défini sur 
un ouvert d'un espace vectoriel de dimension finie ; c'est à cette 
notion de continu (par opposition à « discret », cf. 8 4.2) que notre 
terminologie fait allusion. En réalité, le contexte rendra suffisam- 
ment claire la situation dans laquelle nous nous trouverons, nous 
permettant ainsi de faire l'économie du qualificatif. Nous noterons 
que, dans un cadre topologique très général, les notions de 
« continu » et de « discret » se confondraient. Nous savons que 
l'équation différentielle autonome (E) donne lieu à un système dyna- 
mique. Réciproquement, soit un système dynamique de flot ÿ. Mon- 
trons qu'il est associé à une équation différentielle autonome dont le 
flot est justement Ÿ. 


Raisonnons d'abord par condition nécessaire, en supposant qu'il 
existe f, de classe C K sur un ouvert Q, dont le flot est 4. Nécessaire- 
ment, l’ouvert Q est manifestement égal à l'ensemble des Xo tels 
que (0, Xo) appartienne à %. Soit (Z, X) la plus grande solution au 
problème de Cauchy de condition initiale Xç. On a : 


X(E) = 96 Xo); 


et I est manifestement inclus dans l’ensemble des t tels que (t, Xo) 
appartienne à . De plus : 


1h = 9 
X (©) A TAL Xo)r 
et l'égalité : 
X'(0) = F(X(0)) 
entraîne : 
0 
(0, X0) = f(X0) - 


On voit ainsi que fest nécessairement définie par cette égalité. 
Il y a donc au plus un champ de vecteurs fqui convienne. 


Montrons que fconvient effectivement. Posons pour cela : 
d 
#0) = (0, X5) : 
ot 
cette application est de classe CK sur l’ouvert Q formé des X tels 


que (0, Xo) appartienne à À . Considérons l'équation (E) : X’ = f(X). 
Soit X, dans E. Posons, pour (f, Xç) dans À : 


XD = DCE Xo). 
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On a : X (0) = X. 
De plus : 
1n = 99 
X'() = 55 Xo): 
Or, si (s, Xo) est dans %, en dérivant la relation : 
p CL PES, Xo)) = p(t+ 5, Xo) 


par rapport à t pour t assez proche de 0, ce que permet le caractère 
ouvert de 4, on obtient : 


6 _ 99 
ALLASE X0o)) T st S, X0) . 
Évaluée en t = 0, cette égalité fournit : 
CL) _ 99 
(09 (s, X0)) on TA Xo) , 
soit : 


(E) f(X(s)) = X'(S). 


Ainsi, X est solution du problème de Cauchy ; compte tenu de son 
domaine de définition, elle est clairement maximale ; par consé- 
quent le flot de l'équation (E) est bien égal à 9. 


En résumé, il y a bijection entre les systèmes dynamiques 
continus et les équations autonomes. Il est facile de déterminer 
le champ de vecteurs lorsque le flot est connu. En revanche, la 
détermination du flot à l’aide du seul champ de vecteurs pose 
des problèmes autrement redoutables : c'est à ces problèmes 
qu'est consacrée l'étude des systèmes dynamiques, étant 
entendu que « détermination » ne signifie pas « détermination 
exacte» ni même, dans beaucoup de cas, « détermination 
approchée », mais, le plus souvent, « étude qualitative ». 


4.2 Le cas « discret » 


Nous calquerons la définition d'un système dynamique discret 
sur celle d'un système dynamique continu (8 4.1). Pour cela, nous 
remplaçons la notion d'intervalle de R, qui correspond à un inter- 
valle de temps continu, par ce qui peut tenir lieu d'intervalle de 
temps discret. La plus simple (mais en aucun cas la seule) des pos- 
sibilités offertes est de prendre l’ensemble des entiers naturels. 
Nous sommes donc conduits à introduire la définition suivante. 


Définition 6. 

Soit U un espace topologique. On appelle système dynami- 
que « discret » la donnée d’une application $ de N x U vers U, 
continue par rapport à la seconde variable, vérifiant les proprié- 
tés suivantes : | 


Il 


p (LE(S, Uo)) 
(0, Uo) 


pour tout tet s dans N, et pour tout u, dans U. L'application est 
appelée flot du système dynamique. 


g(t+5, Uo), 


Il 


uo 


Posons : 
f(uo) = O(1, Uo)- 


Une récurrence immédiate montre que : 
®(n, Uo) = F(Uo); 


où ff = fo. o f est l'itérée n fois de fpour la composition. 
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Ainsi, le flot du système dynamique peut être obtenu paritération 
d'une certaine application continue de U vers U. 


Si réciproquement f est une application continue de UÜ vers U, on 
constate que l'on définit un système dynamique discret sur U en 
posant : 


®(n, Uo) = F(Uo). 


En d’autres termes, l'étude des systèmes dynamiques discrets 
équivaut à l'étude de l’itération des applications. 


Une autre façon de présenter les choses est de poser : 
u(n) = f"(Uo), 


définissant ainsi une suite d'éléments de U. On aura alors la 
relation : 


(E) u(n+1) = f(u(n)). 


La suite u satisfait donc une relation de récurrence d'ordre 1, dite 
autonome par opposition à la relation plus générale de la forme 
u(n+1)=gi{n,u{(n)). On peut bien évidemment, et comme dans le 
cas continu, ramener l'étude des suites récurrentes les plus généra- 
les, d'ordre arbitraire, et éventuellement non autonomes, à celles 
des suites précédentes. Nous étudierons désormais l'équation (E), 
lorsque fest une application de U vers U. 


L'analogue du théorème de Cauchy et Lipschitz dans le cas de 
l'équation (E) est le suivant. Soit uw un élément de U. Il existe une 
unique suite u solution de l'équation (E) et telle que u (0) = wo. La 
démonstration de ce résultat découle d’une simple récurrence. 


Nous pouvons à présent, par analogie avec le cas continu, établir 
les principales notions relatives aux systèmes dynamiques discrets. 
On appelle orbite de u, l'image de la suite u (elle-même appelée tra- 
jectoire de uo), c'est-à-dire l'ensemble {u(n)},:N\- Cet ensemble 
est fini ou dénombrable. Lorsque cette orbite est réduite à un 
point &, ce qui signifie que la suite u est constante, ce point vérifie 
f(œ) = &@. On dit que c'est un point fixe de f. Cette valeur sera aussi 
appelée position d'équilibre. 


Si une orbite est finie, il existe p < q tel que : 


u(p) = u(q) = fP-A(u(p)). 


Autrement dit, u (p) est un point fixe de fP"1, et il en résulte que la 
suite (u(n)), sp admet la période p -— q. On dit que la suite u est ulti- 
mement périodique, ou encore périodique à partir d’un certain 
indice. Les valeurs m qui sont une ultime période de u constituent 
un ensemble TN de N\, et Test appelé plus petite ultime période de 
la suite u. La réciproque étant vraie, on voit donc que les trajectoires 
sont ou bien injectives, ou bien ultimement périodiques. 


On pourrait être surpris de ce que l’analogie avec le cas continu ne 
soit pas suffisante pour obtenir, dans le cas d’une orbite finie, une 
suite vraiment périodique. La raison est liée à la forme de l'intervalle 
de temps considéré, qui est N. Si nous remplaçons N par Z, fest 
alors bijective ; les trajectoires finies correspondent bien dans ce cas 
à des trajectoires périodiques. On voit donc apparaître des situa- 
tions particulières liées au type de système dynamique considéré. 
En revenant au cas d'un système dynamique d'intervalle temporel 
égal à N, on remarquera par exemple que deux orbites distinctes ne 
sont pas nécessairement disjointes : on peut seulement affirmer que 
si deux orbites se coupent, leur intersection est encore une orbite. 


Une notion importante est celle d'ensemble stable par le flot. 
Dans le cas discret, un sous-ensemble Vde U'est dit stable par le flot 
lorsque, pour tout n, et pour tout u, dans V, g(n, u,) appartient à V 
Cela revient à dire que, pour tout n: 


FEV, 


ce qui équivaut encore à f(V)E V. 


Cette notion peut d’ailleurs s'appliquer aux systèmes continus. 
Mais, dans ce cas, l’on ne peut en général pas la reformuler à l’aide 


AF 180 - 14 


du seul champ de vecteurs. Il est néanmoins facile de voir, dans les 
deux cas, qu'un ensemble est stable par le flot si, et seulement si, 
c'est une réunion d'orbites. 


5. Quelques exemples 
de systèmes dynamiques 
simples 


5.1 Systèmes dynamiques continus 
dans R 


5.1.1 Étude générale 


Nous considérons ici une équation autonome : 
(E) x = f(x, 


où fest une application de classe C'\ d'un intervalle ouvert Jde R, 
vers R. On notera que le caractère ouvert de J n’est pas essentiel 
dans la suite ; nous ferons cette hypothèse pour nous situer pleine- 
ment dans les conditions de la définition 2. 


Nous appelons Z l’ensemble des positions d'équilibre, c’est-à- 
dire l’ensemble des « de Jtels que f(œ) = 0. Nous limitons notre 
étude aux trajectoires non constantes. 


Soit (Z, X) une telle trajectoire. Puisque ZCJ-Z, I est inclus 
dans une composante connexe /’ de J-Z. Sur J’,f ne s’annule 


pas. Soit F la primitive de L sur J’ qui s’annule en x (0). L’équation 


(E) s'écrit : 
FxXO)X (0 = TT, 


soit F(x(t)) = tt. 
Donc : 
x(t) = Ft) sur I. 


Soit x dans J’. L'application F, définie par : 


est strictement monotone sur l’, et a pour image l'intervalle F(/7’). 
Son application réciproque F7l est donc une application strictement 
monotone sur F(1'), de classe C? sur cet intervalle. 


Si l'on pose x(t) = F-\(t) pour t dans 1’, on obtient donc sur J” 
une solution de (E), telle que : 
X(0) = x, 


qui est la plus grande solution de ce problème de Cauchy. 
Nous avons ainsi déterminé le flot, à l’aide des inverses des primi- 


tives de 5 sur les intervalles où elles sont définies. Nous en dédui- 


sons facilement les orbites. Puisque l'image de F(7’) par F-lest J’, 


on voit que l'orbite d’un élément de J’ est J’. Par conséquent, les 


orbites sont, outre les singletons correspondant aux points d’annu- 
lation de f, les plus grands intervalles sur lesquels fne s’annule pas. 
Nous n'oublierons pas, bien sûr, que la variété des phases est ici 


l'intervalle J, sous-ensemble de R. 
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Dans le cas présent, on remarquera que toutes les trajectoires non 
constantes sont strictement monotones, donc injectives. Cette situa- 
tion très particulière montre que les systèmes dynamiques continus 
de taille 1 ne peuvent guère constituer une source de modélisation 
riche. 


5.1.2 Un modèle biologique 


Pour modéliser l’évolution d'une population x, Malthus propose, 
en 1798, la loi d'évolution : 
x = FX; 
qui conduit à une extinction exponentielle si r< 0, ou à une surpopu- 
lation exponentielle si r> 0. 


Le caractère non réaliste d’un tel modèle conduit Verhulst à pro- 
poser, en 1836, une autre loi d'évolution. Pour cela, il part, comme 
Malthus, de l'existence d’un facteur r, différence entre taux de nata- 
lité et taux de mortalité, qu'il suppose positif. Pour éviter l'inflation 
démographique, il suggère que la loi d'évolution contient un facteur 
de rappel, d'autant plus fort que x est grand. Il est ainsi conduit à 
proposer l'équation différentielle scalaire : 


(F) x = m1 -?) = f(X) 


où fest appelée fonction logistique. 
Les réels positifs Ket r sont des paramètres. Ces paramètres peu- 


vent être éliminés par normalisation. Posant y(f) = kd°) , on voit 


que : 
y = yOC(- y): 

Ainsi, on peut supposer k=r=1. 

L'équation devient : 
(E) x = x(1-X). 

Les équilibres sont «& = 0 et & = 1. Les trajectoires non constan- 
tes sont aisées à expliciter : 

ce’ 


X(T) = : 
Q) 1+cet 


Lorsque x(0) > 0, on obtient le portrait de phase de la figure 8. 


On constate que, lorsque t— +, toutes les trajectoires non 
constantes tendent vers @, qui, pour cette raison, est appelé équili- 
bre asymptotiquement stable. En effet, si la valeur x, est légère- 
ment écartée de la valeur ©, la trajectoire obtenue, qui ne sera plus 
constante, tendra vers la position d'équilibre @. De fait, 
« légèrement » est ici de trop car toute trajectoire tend vers @. On 
exprime cela en disant que le bassin d'attraction de o, est égal à 
R-{@}. 

Quant à o, il n’est limite d'aucune trajectoire autre que l'équilibre 
a : on dit qu'il s’agit d’un équilibre asymptotiquement instable. 


On dispose, dans le cas des équations scalaires, de la possibilité 
supplémentaire de visualiser aisément, non seulement les orbites, 
c'est-à-dire le portrait de phase, mais aussi les trajectoires, en repré- 
sentant le graphe de l’application x. Les équilibres apparaissent 
comme les droites horizontales. La stabilité asymptotique de 1 et 
l'instabilité asymptotique de O0 sont éminemment lisibles sur la 
figure 9. 


Enfin, traçons le graphe de f, lorsque (figure 10) : 
f(x) = x(1- x). 


Les équilibres sont donnés par l'intersection avec l'axe des x. On 
constate que : 


f(a)>0 et f(a)<0. 


Nous généraliserons cette remarque dans le paragraphe 5.1.3. 
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Figure 8 - Portrait de phase dans un modèle avec rappel 
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Figure 9 - Trajectoires d’un modèle biologique dans R 


f(x) 


Figure 10 - Graphe de la fonction logistique f 


5.1.3 Stabilité asymptotique des équilibres 


Nous considérons à nouveau l'équation autonome : 
(E) x = f(x), 


où fest une application de classe C'd'unintervalle ouvertJdeR, 
vers R. Notons Z l'ensemble des points où fs’annule, Z* celui des 
points où fest strictement positive, Z7 celui des points où fest stric- 
tement négative. Soit & un élément de Z. 


AF 180 -15 


PRÉSENTATION DES SYSTÈMES DYNAMIQUES 


Proposition 4. 
Dans les conditions précédentes, on a les assertions suivantes. 


æ Si f(a)<0, a est un équilibre asymptotiquement stable. Soient 
la’, al et la, &”’[ les composantes connexes de J - Z dont æ est 
une extrémité. Le bassin d'attraction de « est alors ]a’, &’”T. 


m Si f(a)>0, « est un équilibre asymptotiquement instable. Son 
bassin d'attraction est réduit à {œ}. 


æ Si f(œ) = 0, æ est un équilibre asymptotiquement stable si, et 
seulement si, à est extrémité gauche d’une composante connexe de 
J-Z* et extrémité droite d’une composante connexe de J-ZT. Le 
bassin d'’attraction de a est alors la réunion de ces deux composan- 
‘ tes connexes et de {æ}. 


Preuve. © Nous prouverons le premier point, les deux autres étant 
très analogues. . 


Si f(æ)<0, il existe un voisinage à gauche épointé de & sur 
lequel fest strictement positive. Soit ]æ”, al la composante connexe 
de J-Z contenant ce voisinage. Reprenons les notations du 
paragraphe 5.1.1 ; l'intégrale qui définit F est divergente en a. Par 
conséquent, F (]æ’, al) contient [0, +[ et Ftend vers + en «œ. 
Donc F1, c'est-à-dire la trajectoire x, est définie sur [0, +[ et 
tend vers œ7en +. 

Un raisonnement analogue à droite montre que la’, "TI est 
inclus dans le bassin d'attraction de a. 


D'un autre côté, les orbites qui ne sont pas l’un des deux interval- 
les la’, al ou la, æ’’[ sont à une distance strictement positive de &. 
Les trajectoires correspondantes ne peuvent donc tendre vers & © 


5.1.4 Exemple de système dynamique 
avec paramètres 


L'étude de la prédation d'un épicéa par des larves conduit, après 
normalisation, à proposer comme modèle d'évolution de la popula- 
tion d'arbres l'équation différentielle autonome : 


. 2 
(E) x = 1 7) = £.x00. 


Les réels strictement positifs r et k sont des paramètres qui vont 
ici nous intéresser au premier chef. La normalisation a permis de 
réduire le nombre initial de paramètres de 4 à 2, sans toutefois les 
supprimer complètement. 


5.2 Systèmes dynamiques discrets 
dans R 


5.2.1 Généralités 


Soit fune application continue d'un intervalle Jde R,, vers J. Pour 
u, dans J, on considère la suite u telle que : 


ü(n) = (0); 


Nous nous intéresserons ici particulièrement au comportement 
asymptotique de cette suite lorsque ntend vers +.Remarquons, 
d'abord, que si la suite u converge dans J, sa limite & ne peut être 
qu'une position d'équilibre, par passage à la limite dans l'égalité : 


u(n+1) = f(u(n)). 


Un équilibre symptotiquement stable (appelé encore point fixe 
attractif) ou instable est défini comme dans le cas des systèmes 


Figure 11 - Visualisation de l'orbite d'un point 


dynamiques continus. L'ensemble des u, tels que (u (n)) tende vers 
æ est appelé bassin d'attraction de a. 


La recherche de points fixes est souvent un préalable à l'étude 
d’un système dynamique discret. Notons à ce sujet le résultat élé- 
mentaire suivant. 


Proposition 5. 


Soit f une application continue de source le segment J, et à 
valeurs réelles. 


& Si Jcf(J), f admet au moins un point fixe dans J. 


& Si J>f(J), f admet au moins un point fixe dans J. 
Preuve. © Posons J= [a, b], avec a < b et @(x) = f(x)- x. 


Dans le premier cas : 
(a) <0 et #(b)=0, 


donc ÿ s’annule entre a et b. 
Dans le second cas, soient ce J tel que f(c)=aet de J tel que 
f(d)=b.Ona: 
(c) = a-c<0 et ÿ(d)= b-dxz0, 


donc ÿ s'annule entre cet d. 
© 


Pour étudier un système dynamique discret dans R,, il est souvent 
commode de faire figurer le graphe de f, ainsi que la première bis- 
sectrice. Les points d'intersection correspondent aux équilibres. Par- 
tant d’une valeur u, on peut placer les valeurs u (n) en remarquant 
que u (n + 1) est obtenu en projetant le point P, sur l'axe des x. 


La figure 11 correspond à f(x) = -sinx et us = 0, 5. 


NI 
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5.2.2 Cas des applications strictement monotones 


Nous étudions dans ce paragraphe un cas particulièrement sim- 
ple, celui où l'application f, toujours supposée continue, est une 
bijection d’un intervalle J de R sur lui-même ou, ce qui revient au 
même, une application strictement monotone. En réalité, cette situa- 
tion correspond à un système dynamique discret dont l'intervalle 
temporel est Z (cf. 8 4.2). 


Pour la commodité, nous étudierons le cas où J'est ouvert, le cas 
général étant obtenu par de légères modifications des raisonne- 
ments. 

# Cas d’une application strictement croissante 


Nous notons Z l'ensemble des points fixes de f; soient ug un élé- 
ment de J -Zet I la composante connexe de J - Zcontenant ug. On 
constate que Z est stable par f. En effet, soit xe I ; notons 
a = sup. 


Si &œ = +, on a bien sûr f(x) < à. 
-Si œ est un point fixe de f, x < œ implique : 


F0) <a) = «à 


car fest croissante. 
Enfin, si æ est un réel qui n’est pas un point fixe de f, c'est que 
a = sup J, et comme f(J)cJ, on a encore: 
f(x) < à où F(X)< a 


selon que «æ est dans J ou non. 
On raisonne de même avec le borne inférieure de 1. 
Sur 1, f(x) -xne s'annule pas, donc on a toujours : 


f(xX)> x ou F(X)<xX. 


Dans le premier cas, on voit par récurrence que la suite u est stric- 
tement croissante et, dans le second, qu'elle est strictement décrois- 
sante. 


Plaçons-nous dans le premier cas. Si la suite est convergente dans 
J, sa limite «& est un point fixe de f, qui ne peut être que la borne 
supérieure de Z. Si réciproquement sup J est un point fixe de f,uest 
majorée, donc converge, vers Sup 1. 


En résumé, lorsque f(u,) > u,, la suite u est strictement crois- 
sante : elle converge dans J si, et seulement si, u’ est majoré par 
un point fixe de f, et sa limite est alors le plus petit point fixe de 

| fsupérieur à üo. 

| On obtient un résultat analogue dans le cas où f(u,)<u,,en 

| renversant l’ordre. 


La figure 12 correspond à l'application f définie par : 


f(x) = V/x+2 


strictement croissante sur R,. Cette application admet l'unique 
point fixe 2. Pour ug < 2, la suite (u (n )) converge donc vers 2. 


Cas d'une application strictement décroissante 


L'application g = fo f est alors strictement croissante. On peut 
donc étudier séparément les suites (u (2n)) et (u (2n + 1)), solutions 
de la récurrence : 


v(n+1) = g(v(n)) 


avec, respectivement, les conditions initiales v(0)=u et 
v (0) = f(Uo). 


0 0,5 1,0 1,5 2,0 2,5 


Figure 12 - Visualisation d’une orbite (f croissante) 


Pour terminer l'étude, il suffit de remarquer que la suite u con- 
verge dans Jsi, et seulement si, les suites (u (2n)) et (u (2n + 1)) con- 
vergent dans Jet ont même limite. 


La figure 13 correspond à l'application définie par : 


F9 = 34-22), 


. strictement décroissante sur R, . Cette application admet l'unique 


point fixe 1. Pour up = 0, les suites (u (2n)) et (u (2n + 1)) convergent 
toutes deux vers 1. Donc (u(n)) tend vers 1. 


Figure 13 - Visualisation d'une orbite (f décroissante) 
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